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practicality of replacing costly general purpose dieactal 
computers with relatively inexpensive dedicated micro- 
computers on-board naval aircraft. The microcomputer 
program showed that microcomputers have sufficient speed 
and accuracy to solve the navigation Problem. inorder 
to overcome the microcomputer's major deficiencies, speed 
and accuracy, special arithmetic subprograms based on 
table look-up were developed to trade inexpensive memory 
for more speed. An application of graph theory in the form 
of process graphs was made to facilitate the development 
and documentation of the navigation program. To aid. in 
the testing of the microcomputer program, a Fortran simu- 
lation program was developed to confirm the results of an 
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ieee LNT RORUCTION 


Current naval aircraft are depending more and more on 
airborne digital computers. The digital computer has proven 
to be a great aid, providing the precision and speed required 
to perform many calculations. The digital computer is very 
versatile, because the same computer can be programmed to 
perform new or different tasks. The only limitations to 
the use of digital computers are those of cost and 
maintenance requirements. 

The digital computers currently used by such naval 
aircarft as the P3C, A6E, ATE, and S3A are all large general 
Purpose computcrs. These computers are very fast, flexible, 
but expensive. The large size, large power requirements, 
and great cost has limited each aircraft to one Such com- 
puter. The complexity of these computers has made main- 
tenance difficult. The large cost of each unit makes spares 
prohibitive, therefore one computer going down results in 
wiewoperatvional oss of one aircrait . 

It is advantageous to utilize several small distributed 
computers to meet all of the systems requirements. In 
order to minimize the probability of the entire system 
failing due to the failure of a single component in a cri- 
Mica MecOnputer Cireurl, complevely separate computers could 
be used for the various system requirements with back up 


computers ready to fill in when needeea. In addition to 


14 





reliability improvements, a distributed approach offers 
icmeessltbliatyeet USing less complex equipment, flexibility 
in matching equipment to system requirements, and increased 
Standardization. The limitation of this approach is the 
increased cost, weight, size and the complexity of system 
Picemeconniect ons. 

The -creation-of -the-mtcrocomputer, using new developments 
in the Large Scale Integration (LSI) technology, has made 
the distributed computer system possible. The microcomputer 
is a complete general purpose computer on a set of four 
Standard LSI chips. The LSI chip measures 200 mils by 200 
mils, requires less than one watt of power and costs about 
$30. The processors Central Processing Unit (CPU), constructed 
On a Single chip, is designed to be used in a multiprocessor 
eutronmens. —~ lhe limitations of a2 microcomputer are 4 
limited instruction set and slow speed. 

Mis thesis 1s°a design study of the possibility of 
using the MCS-4 microcomputer as the Avionics Navigation 
Computer in a complex navigation system. Section II dis- 
cusses what a microcomputer is, what microcomputers are 
currently available, and what is the system makeup of the 
MCS-4 microcomputer used in this report. Section III des- 
eribes the navigatin systems currently available, what their 
advantages and limitations are, and the navigation system 
Gm@een for this report. The navigation equations used to 
compute the current position of the system's carrier are 


discussed in Section JV. The actual propramming of the 
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Plcmelecnmolvciy ts discussed in deétall in Section V. This 
dneludes the programming ailds developed in this report, a 
program analysis of the requirements of this system, a 
detailed discussion of the executive routine and the sub- 
routines used in this program, and an error bound analysis 

of the final program. Section VI describes the FORTRAN 
‘Simulation “program written’to aid in the writing, debugging, 
and testing of the final navigation program. The conclusions 
of this report are summed up in Section VII. The listing 

of both the MCS-4 microcomputer navigation program and the 


FORTRAN simulation program are included in the Appendices. 
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iieiLCMOCOMEUTER 


A microcomputer is a general-purpose digital computer 
Poncumcted trom a set of LSI chips. It has a complete 
instruction set and iS capable of addressing sizeable memo- 
ries. It can interface with a full complement of input and 
output devices. 

The main SOAS of a microcomputer is the micro- 
processor. The microprocessor is a CPU on a chip that 
interprets and executes instructions in a bit-parallel 
fashion. Included on the CPU chip are the Index registers, 
| Arithmetic unit and Input/Output control logic. 

fiaewieain teacure thay diGtancuishes 2 miercorocessor 
from a general purpose or minicomputer CPU is that the 
entire CPU is on one chip. This has been made possibly by 
Large Scale Integration, where over 14,000 Metal-Oxide Semi- 
eonductor transistors can be put on one chip. A minicomputer 
CPU with the capability of one microprocessor would require 
evere OO TTL packages. The major advantages of micropro- 
eessors are low cost, low power requirement, and less com- 
plexity in system design. Figure 1 is a photomicrograph 
of the 4004 CPU chip used in the MCS-4 Microcomputer. 

The microprocessor becomes a microcomputer when a con- 
trol program, memory, and input/output circuits are added 
CO the system. The control program 1S usually Metal—Masked 


on a Read-Only Memory (ROM) chip, however, a Programmable 
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Read-Only Memory (PROM) may be used for initial system 
development. The storing of data or variables in a program 
4s handled by the Random Access Memory (RAM) chip which is 

a read/write type memory. The input/output ports are incor- 
porated on the RAM, ROM, or CPU chip or any combination 

of these chips. The total chip area of a Cyc cles Iise 2 O— 
computer is 200 mils by 200 mils and costs less than $200. 

Microcomputers have been used in specialized, single 
user data processing systems and as components in digital 
products. Microcomputers are currently being used in cal- 
culators, special purpose terminals, measurement systems, 
intelligent traffic controls, small business computers, and 
digital cash registers. Research is in progress to use 
microcomputers in intelligent peripherals, multap lexors 
and communications controllers, automotive control systems, 
and educational systems. Microcomputers can also be used 
in large computer systems to relieve the large central 
processor of much of the overhead associated with scheduling, 
text editing, and file management. 

Microcomputers major application will be in the area of 
dedicated computations. They provide all the Gata processing 
power needed for these apolucaceons.w LHese areas include 
computational tasks required by NTDS, TSC, and most current 
mavalfeailrvcrat - 

The following is a list of advantages of MLerocompucers.: 


1. Reduced costs due to reduction in number of LGELe 
Camdarey DCS. 
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. Self test capability. 
. Equipment modularity. 


. Equipment commonality. 


2 
3 
4 
5. Ease in design changes through microprogram. 

6. Reduced logistic support. 

if Standardization of peripheral interfaces. 

8. Multiprocessor capability. 

The total savings in cost, power required, and size, 
together with the flexibility inherent in the systen, make 


the microcomputer a powerful tool for system design of 


dedicated computational tasks. 


A. COMPARISON OF MICROCOMPUTER WITH GENERAL PURPOSE COMPUTER 

The general purpose computer and microcomputer represent 
the two extremes on the computer scale. The general purpose 
computer is large in size, very expensive, requires much 
power, and is very fast. The microcomputer is small, inex- 
pensive, requires little power, and is relatively slow in 
comparison. The general purpose computer is built using 
Integrated Circuits (I.C. Technology), while the microcom- 
puter is built from Large Scale Integration Techniques Clon 
Technology). The great cost savings of LSI Technology in 
computer design can best be shown by comparing the cost of 
designing and manufacturing a 3000 gate logic Clue ewes 1s 
I.C. techniques versus the LSI techniques. 

The cost of ManuTpacluring a computer system can be 


expressed by the following formula: 


20 





MCChCT) teers) VC) + CGC) | 
al 


© 

tA 

ll 
lms 


where 


CS = Cost of System 

CM(i) = Cost of Manufacturing the i-th Replaceable 
Module 

CT(i) = Cost of Testing the i-th Replaceable 
Module 

V(i) = Number of Modules used in the system of 


the i-th type 


CG(1) = Cost of generating the i-th Replaceable 
Module, including design, layout, test 
condition generation, etc. 


The cost of manufacturing the computer system from 1.C. 
circuits was estimated using figures from people knowledgeable 
in this area, however, documentation of these figures was 
not available. 

The cost of manufacturing and the cost of testing an 
individual I.C. module was assumed to be independent of the 
type of I.C. module because of standard means of manufac- 
turing and testing these units. This assumption simplifies 


the "cost of system" equation as follows: 


Me 


CMCi) CM 


Me 
c 
= 


Cea) 


eal: 





where 


V(i). 
at 


<j 
0 
wus 


a 


Z.C, technology does not permit highly complex functional 
circuitry to be placed on one module. The typical module 
consisting of sixteen pins contains elght gates. A 3000 
gate logic pattern would therefore require V = 360 modules. 
The technology typically used by computer manufacturers 
permits 60 modules to be assembled on one panel which pro- 
vides all the interconnections between the modules. The 
system would then require six panels at an estimated cost 
of $100 per panel. 

In determining the number of distinct modules used, it 
was‘ assumed that about one-half of the modules used in the 
system are duplicates of previously designed modules, n = GOr 
The cost of designing one new module was estimated to be 


$250. The total cost of one system was computed as follows: 


CM = $20 
CT = $5 
V = 360 


CG(i) = $250 
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n = 180 


CS = (CM + CT) + CG(i) + 600 


CS 9,000 + 45,000 + 600 = 54,600 

To produce this unit in large numbers, the cost of manu- 
‘facturing and testing the modules and panels is multiplied 
by the number of units produced while the cost of designing 
remains unchanged. The cost of 1000 units could be predicted 


as follows: 


CS(1000) = 9000(1000) + 45,000 + 600(1000) 


Tora cost of 1000 merits =39-645 5000 


Next the cost of developing and manufacturing a 3000 
gate logic unit using LSI technology was investigated. The 
cost estimates for this analysis were obtained from a trip 
PoOminte  iCOrpa. oanta Clara. California... The Lol design 
cost of the 8008 CPU chip used in the MCS-8 Microcomputer 
was used as an indication of the cost associated with manu- 
Pevunin es bne reduared 1OR1¢C nite by Lol technology. 

The entire design and layout of the 8008 CPU was done 
Sr,elusavely by hand and took a total of four man-years. 

The approximate cost of developing, this chip is estimated 


ac 100.000, “The total manufacturing and testing costs of 


2 3 





the chip is $15. Since all the logic gates are on one chip, 
there are no other costs associated with the Unite » che 


cost of one system is therefore: 


CS Cia) ol) + LOO. OCG 4 0 


cS = $100,015 

Again, the cost of producing the system in large quan- 
tities, the design cost remains fixed while the manufacturing 
and testing costs are multiplied by the number of units. 


The cost of 1000 units using LSI technology is as follows: 


CoGoCe e— lS) ClUGO ee A00 00u 


Total cost of 1000 units = $115,000 


It was concluded from these calculations that LSI 
technology could produce the system in quantity about 80 
Gfimes cheaper than by I.C. technology. 

It was concluded from this analysis that the microcom- 
puter has a tremendous cost advantage over peneral purpose 
computers in systems produced in large quantities with a 
need for a dedicated computer. The object of this thesis 
4s to demonstrate that the microcomputer also has the speed 


and computational power to handle complex tasks. 
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B. COMPARISON OF MICROCOMPUTER WITH MINICOMPUTER 

The minicomputer was developed to meet the need for 
dedicated systems to handle data processing requirements 
such as communications control, data acquisition, and small 


business accounting. The microcomputer is proving itself 


Be — — 


capable Onl handling these tasks with improved price/perfor- 
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mance , compactness, and reliability. 


The greatest advantage the minicomputer has over the 
microcomputer is speed. The LSI chips used for microcom- 
puters are made using Metal Oxide Semiconductor, MOS, tech- 
nology. MOS technology allows smaller size of individual 
transistors, and logical circuits together with low power 


consumption. The electrical properties of MOS circuitry, 


a pe 
i —— — a 


however, make it slower than the Bipolar circuitry used am 
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minicomputers. In order to increase the parce ecionnl speed 
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ster as ge oT ee ole the computational “speed of microcomputers 
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can become as fast or caver then Se 


= ——— a —— 


The Ree. aeeare the ee ox: has over the micro- 
computer is in the instruction set size. The first genera- 
tion microcomputers were limited to between NO to 60 instruc 
tions while the minicomputers had instruction sets in the 
100-120 Sead The gap here is also closing with the second 


generation microcomputers having instruction sets of Oe U0 


instructions. 
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The last major advantage of the minicomputer is the 


existing software developed over the past several years to 


assist program development. h the development of high 


level programming languages (PL/M for the Intel 8008 ana 
ee ee a eee ae 
8080), the microcomputers are quickly eliminating this 
Se 
advantage of the minicomputers. 
Ls  t—“<OCPr”rt”tttrrmrmlmlml 


oe 





The microcomputer offers better price/performance, 
lower power consumption, smaller size and higher reliabil- 
ity than minicomputers. Although a Single microcomputer 
can not match the power of a minicomputer, Several micro- 
computers can be combined to share the workload at a cost 
still less than a minicomputer. An added advantage of 
multi-microcomouters is processor reliability which can be 
increased through the use of back-up processors, providing 
a self-test and repair capability. 

Table I is a summary of the comparison of microcomputers 


and minicomputers. 


C. CURRENTLY AVAILABLE MICROCOMPUTERS 

The number of microcomputers being designed and currently 
available is increasing rapidly. This section covers eight 
microcomputers that are currently available or will be in 
Bae near future. Tabdie Il is a Summation of the capabilities 


and support available for each microcomputer. 
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CPU Cost 
Instruction Speed 


Execution Speed 
(Memory to Memory Add) 


Instruction set 
“Price/Instruction Ratio 
Registers 
Price/Register Ratio 
Memory Size 


Performance/Price 
Word Length 


Add Time x Price 


Power Consumption 


meliability (CPU) 


MINICOMPUTER 


$1500 
2-5 msec 


5-20 msec 


100-150 
$80-$300 
1~30 
$300 

64K 


200 (Bits 
Sec-pol lar) 





4 Watts 


Less (Due to 


100 TTL Packages) 


MICROCOMPUTER 


$30 
10-20 msec 


15-60 msec 


50-100 
$40 
16-100 
$15 
64K 


hOOm CB acs 
Sec-Dollar) 





1 Watt 


Greater (Due to 
Gnie Lol ens 


TABLE I. MICROCOMPUTER VERSUS MINICOMPUTER 
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ills Intel MCS-4 

The MCS-4 is a four-bit parallel processor with a 
fixed instruction set. The 4004 CPU is a P-Channel MOS LSI 
chip mounted on a 16-pin package. A four-bit data bus 
connects the CPU with up to 16 ROMs and 16 RAMS. The 
instruction cycle is 10.8 microseconds. 

The instruction set used by the MCS-4 consits of 
4S instructions grouped into three categories: Machine, 
Input/output, and accumulator. System development aids 
include a cross assembler available on the IBM 360/67 located 
at the Naval Postgraduate School. The SIMQ4-02 Hardware 
Prototyping System with Assembling, Programmable-ROM 
Programming, and Debugging capability is also available 
ee che NPS microcomputer lab. 

The MCS-4 Microcomputer was chosen Necene navigation 
microcomputer in this thesis. The MCS-4 is covered in more 
depth in Section II.D. 

2. Intel MCS-8 

The MCS-8 uses the 8008 single CPU chip. The 8008 
fea elpht—bit fixed instruction set parallel processor 
mounted in an 18-pin package. The 8008 executes a single 
instruction in 20 microseconds. The 8008 is capable of 
addressing 16K bytes of memory. An eight-bit data bus 
interfaces the processor with memory. A total of 48 instruc- 


oilers are avaitabple broken into four groups: Instruction 





Register, Accumulator, Program Counter and stack control, 
and Input/Output. 

System Development aids include a cross assembler 
available on the IBM 360/67 at the Naval Postgraduate School. 
The SIM8-01 Prototyping System with programmable and erasable 
ROMs is available in the NPS microcomputer lab for develop- 
ment and check-out of microprograms. The development of 

the PL/M Higher Level Language patterned after PL/J for 
the MCS-8 greatly facilitates the programming task. the 
PL/M compiler is also available on the IBM 360/67 at NPS. 
Lastly the INTELLEC-8 system, available in the NPS micro- 
eomputer lab, makes available to the programmer a resident 
software monitor, assembler, PROM programmer, and text 
e@auing capability to aid in the Gevelopment of micropro- 
grams for the MCS-8. 

3. Intel 8080 ee 

The Intel 8080 is the first of the second generation 
microcomputers. The 8080 derives more speed and capability 
by using the more efficient N-Channel MOS. The 8080 CPU 
4s a single-chip eight-bit parallel processor in a HOmpan 
package. The 8080 contains six eight-bit data registers, 
an eight-bit accumulator, three eight-bit temporary regis- 
ters, four testable flags, and an eight-bit arithmetic/logic | 
unit. The execution time for one instruction is 2 micro- | 
seconds. 

The 8080 can directly uccess 64K bytes of memory. 


A separate 16-bit address bus is provided as well] as ten 
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control lines that indicate CPU and I/O bus status. Up to 
256 I/O devices can be directly addressed. Multiprocessor 
capability is designed into the MCS-8080. 

The 8080 is software compatible with the 8008 micro- 
processor. The 8080 instruction set contains the 48 instruc- 
Euems Of the 8008 plus 26 new instructions for a total of 74. 
The 8080 is capable of unlimited subroutine nesting. 

system development aids include a cross assembler, 
INTELLEC 8080 simulator, and the PL/M higher level program- 
ming language. 

4, AMI 7300 ee 

The AMI 7300 is an eight-bit fully parallel, bus- 
oriented processor. The processor consists of two chips, 
the Micro instruction ROM chip (MIR) and the registers- 
adder logic unit chip (RALU). Both chips are P-Channel MOS 
and each is packaged in a 40-pin dual-inline package (DIP). 

The MIR contains a mask-programmable 512 word x 
fe bit microinstruction ROM and a programmable instruction 
mapping array, allowing up to 50 microprogram locations to 
mes predefined for macroinstruction decoding. This allows 
momiasS>ruction set to be tailored to suit a particular 
application. A hardware address stack and loop counter 
allows subroutine nesting to seven levels. 

The RALU contains 48 registers which may be utilized 
a@owone or two first-—in/last-out stacks or as one or two 
files of general registers. The eight-bit adder/subtractor 


performs over j0 arithmetic and lorical operations. The 


of 





instruction set consists of three basic formats: Register 
control, literal and branching. The processor can address 
up to 64 K of memory. 

Software development aids include a cross assembler 
and instruction simulator. 

5. Fairchild PPS-25 

The -PPS—-25 is a BCD oriented 25-digit .-serial/paralilel 
processor. The system is best suited to decimal applications 
such as calculators, keyboard/printer interface, and vending 
machines. Four level subroutine nesting and three way condi- 
minal branching are provided. The system can be micro- 
programmed with a custom instruction set to best meet the 
needs of the user. 

Program storage consists of up to 26 ROMs, each 
capable of storing 256 twelve-bit words. Seven general 
purpose 25-digit registers are provided and an external 
iieerrupt capability is inciuded. 

The 3805 Arithmetic Chip includes the Adder/Subtractor 
plus a 25-digit register. The instructions are located in 
the 3810 ROM. A total of 30 arithmetic/logic and Out © 
instructions are available with the standard set. The I/0 
format permits expansion to 63 I/O instructions. 

Software development aids include a cross assembler 
aad instruction simulator. 

6. National Semiconductor GPCP/IMP-16 
The GPCP sixteen-bit microcomputer processor consists 


of five MOS LSI chips, each mounted in a 40-pin DIP. The 
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five chips consist of four RALU chips and one CROM (Control 
Read Only Memory). Each RALU chip is a four-bit slice of 
CPU with its own registers, ALU logic, and I/O data lines. 
Multiple CROMs may be used to increase the Size of the 
Migatruction set. 

The RALU consists of a four-bit program counter, 
four-bit memory data register, four-bit memory address 
register, four-bit accumulator, a pushdown stack, data 
multiplexer, and four-bit arithmetic and logic unit. The 
system may be expanded to 32 bits by adding four more RALU 
ehips. 

Hine weno eecentains the Control ans@erucctons tor the 
RALU chips. The CROM is broken into two parts, an instruc- 
tion ROM containing 100 twenty-three-bit words and an address 
control ROM consisting of 12 programmable ten-bit words. 

The IMP-16 is a sixteen-bit microcomputer developed 
momuse the GPCP microprocessor. The CROM provided with the 
system contain a 43 word instruction set. The instruction 
set may be expanded to meet the system designer's specific 
needs. Communication between the RALU chips and the CROM 
chip is over a sixteen-bit data bus and a four-bit control 
bus. This requires sixteen extra TTL packages excluding 
memory and timing. 

Software development aids include a cross assembler 
ara prototyping: system with resident monitor, assemblers, 


and linking loader. 
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7. North American Rockwell PPS-4 

The PPS-4 microcomputer consists of a set of six 
MOS Lsk chips. The CPU is a four-bit single chip processor 
mounted in a 42-pin package. The other five chips support 
the CPU and consist of: A 256 x 4-bit RAM, 1024 x 8-—bit 
RAM, an I/O buffer, and a two-phase clock generator. 

The CPU can drive up to 4K bytes of ROM and 4k 
bytes of RAM over its 12-bit parallel re bus. The 
Pas Hem INnShrcuchlon Sel COntains 50 instructions with an 
execution time of five microseconds. 

Twenty-one multiplexed lines interconnect thes cru 
with ROM, RAM, and I/O circuits. These lines are func- 
tionally grouped into twelve parallel address lines, eight 
parallel data lines, and one write command and I/0 enable 
line. The ROM has two chip-select inputs and the RAM has 
one chip-select input, which may be directly controlled by 
Mesechere Ourpurs from the CPU or 1/0 circuits to expand on 
memory without the need for auxiliary circuits. Each I/0 
fio canenandlevup, tO 12 inputs and 12 outputs with a total 
Syatem Capability of 16 2/0 chips. 

Software development aids include a cross assembler 
and Simulator available on a national time-sharing network. 

S. Sigpetics PIP 

The programmable integrated processor (PIP) is a 

single-chip eight-bit CPU made with N-channel MOS technology. 


This second generation CPU is packaged in a 40-pin DIP. The 
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address logic, control memory, and ALU are organized around 
paeeteht—bit bidirectional data bus. There are fifteen 
address lines for handling external memory and I/O circuitry. 
The PIP instruction set contains over 64 instructions with 
an execute time of less than ten microseconds for the most 
eemolex instruction. 

Wao wir Chtpe carpe” oroken™into"four parts, the 
podress logic, the RALU, and the control section. Ihe 
address logic section handles all instructions and includes 
a return address stack that allows subroutine nesting to 
eight levels. The RALU section contains four 8-bit general 
purpose registers and executes all arithmetic, boolean, 
compare, and rotate operations. The control section manages 
operation of all external control lines, decodes all instruc- 
mies, and cOOrdinates the activities of all other internal 
etreultry. 

Software development aids include a cross assembler 
pm@eanstruction simulator. 

9. Summary 

Microcomputers have the capability of replacing 
both special function logic modules and large computational 
machines. Certain microcomputers are more suited for one 
moolication than for others. In selecting a microcomputer, 
Such parameters as data word length and type, instruction 
power, and interface structure must be considered. Table Il 


is a summary of the microcomputers discussed in the section. 





The rapid developments in the LSI technology should 
bring about many new improvements in microcomputers such as 
wider word lengths, larger memory capacity, and more flexi- 
ble and convenient I/0 interfacing. Improvements in LSI 
structure will make possible single chip microcomputers with 
the CPU, I/0, and memory all on one chip. Work on a Bipolar 
LSI processor could develop into a microcomputer with more 


Speed and memory than most present minicomputers. 


D. THE MCS-4 MICROCOMPUTER 

The MCS-4 was the first microcomputer made. Compared 
memocher Microcomputers, it 1s Slow with a limited instruc- 
tion set. The most complex function available on the MCS-4 
is a four-bit add. The MCS-4 is also the least expensive 
microcomputer. The structure of the MCS-4 is similar to a 
general purpose computer, making it compatible with the 
requirements of a navigation computer. 

The MCS-4 was chosen as the microcomputer in this design 
Study for two major reasons. The first reason is that the 
McS-4 is the least powerful and hence serves as a lower 
bound of the microcomputers. To prove that the MCS-4 is 
capable of handling the required navigational computations, 
would in itself prove microcomputers capable of handling 
complex tasks. The Second major reason for choosing the 
MCS-4 microcomputer is that it is available, has been tested, 


and has the required software aids to complete a design study. 
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1. McS-4 System Description 

The MCS-4 microcomputer is built up from a standard 
Seteot off—-the=shelf chips. The only custom part is the 
ROM chip which stores the specific program defined by the 
user and requires a metal mask option for each new program. 

The MCS-4 consists of four chip types, each packaged 
ima conventional 16—-pin DIP: 

(1.) A Central Processor Unit Chip-CPU-4004 

(2.) A Read Only Memory Chip-ROM-4001 

(3.) A Random Access Memory Chip-RAM-4002 

(4.) <A Shift Register Chip-SR-4003 

The CPU contains the control unit and the arithmetic 
unit. The ROM stores the program and data tables, the RAM 
stores input data and variables, and the Shift Register is 
used in conjunction with I/0 devices to effectively increase 
the number of I/O lines. 

A complete microcomputer can be built using only 
a single CPU chip and a single ROM, and the only external 
eeroultry requirec is a two phase clock. The CFU is capable 
of driving a system up to 16 ROMs (4K bytes), 16 RAMs (640 
Pees) ana 128 1/0 lines, with no additional interfacing 
circuitry. The CPU communicates with the RAM's and ROM's 
iy means of a four-line data bus. This single data bus is 
used for all information flow between the chips except for 
control sigznals which are sent to RAM and ROM over five 
additional lines. Figure shows the MCS-4 System Inter- 


connection. 
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The MCS~4 ee. a 10.8 microsecond instruction cycle. 
The cycle is broken into eight steps. In the first three 
steps, the CPU sends the memory address to the ROM in three 
4-pit nibbles. The ROM then sends back 8 bits of instruction 
in two 4-pit nibbles during steps 4 and 5. The instruction 
is then interpreted and executed by the CPU during the last 
““Phypee “steps. 
a. 4004 CPU 
The heart of the MCS-4 microcomputer is the 4004 
CPU. The 4004 CPU contains the following functional blocks: 
(1.) Address Register and Address Incrementer 
(2.) Index Register 
(3.) 4-pit adder 
(4.) Instruction Register, Decoder and Control 
(5.) Peripheral Circuitry 
The Address Register is a RAM array of 4 x 12 
bits. One level is used to store the current instruction 
address, leaving three levels to store the addresses of 
nested subroutines. AS each byte of address is sent onto 
the data bus, the address is incremented by a 4-bit carry 
look-ahead circuit. The incremented address is then trans- 
ferred back into the address register. 
The Index Register is a RAM array of 16 x 4 
bits and has two modes of operation. In one mode of opera- 
tion the index register provides sixteen directly addressable 


storage locations. In the second mode, the index registers 





provide eight pairs of addressable storage locations for 
addressing RAM and ROM as well as for storing data fetched 
from ROM. The index registers can thus provide 64 bits of 
RAM to a minimum MCS-4 system of one CPU and one ROM. 

The H-pit Adder is of the ripple-through carry 
type. The output of the adder is transferred to the accumu- 
lator and carry flip-flop. The accumulator is provided with 
a shifter to implement rotate right and rotate left instruc- 
tions. The accumulator also communicates with the command 
control register, the condition flip-flop and the 4-bit 
internal data bus. The condition logic allows the execution 
BGeeconditional instructions based on the @eontents of the 
accumulator, index registers, or the status of the control 
lines. 

The Instruction Register is an 8-bit register 
which is loaded with the two 4-bit nibbles of instruction 
mead from the ROM. The instructions are decoded in the 
instruction decoder and appropriately gated with timing 
signals to provide the control signals for the various 
munetLOonal blocks. 

The peripheral circuitry consists of the 4-bit 
internal deta bus, the timing and SYNC generator, one ROM 
command control and the four RAM command control output 


Burfers, and the resét flip-flop. 
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b. 4001 ROM 
The 4001 is a 2048-bit metal mask programmable 
ROM. The 4001 performs two functions. As a ROM, its first 
function is to store 256 x 8-bit words of program or data 
tables. The second function of the 4001 is to act as a 
vehicle of communication between the data bus lines and 
peripheral devices through the LSpeee eel tb —-OutouL DOLL 
located on each chip. 
ec. 4002 RAM 
The 4002 also performs two functions. As a 
RAM it stores 320 bits arranged as four registers of twenty 
h-pit characters each. As a vehicle of communication with 
oeripheral devices, it is orovided with LM outowT lanes cald 
associated control logic to perform output operations. 
d. 4003-SR 
The 4003 is a 10-bit serial-in/parallel-out 
serial-out shift register. Its function is to increase the 
number of output Jines to interface with I/O devices. 
2 EinstuNer lon, Sel 
The MCS-4 Instruction Set consists of a total of 
45 instructions grouped into three sets: Machine iInstruc- 
tions, Input/Output and RAM Instructions, and Accumulator 
Group Instructions. 
The Machine Instructions are the Housekeeping 
fmstructions Of the IMCS-}. They consist of two types of 


instructions: one-word instructions which are Spine wilde 
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and require one instruction cycle, and two-word instructions 
whieh are 16-bits wide and require two instruction cycles. 
Table III is a list of the Machine Instructions. 

The Input/Output Instructions are used to transfer 
data between the Accumulator and RAM. These instructions 
are also used to transfer data between the Accumulator and 
the I/O ports located on the ROMs and RAMs. Table IV is 
meeiccect the Lnput/Output Instructions. 

The Accumulator Instructions are used to perform 
mit=by-bit manipulation of the data in the Accumulator. 
Table V is a list of the Accumulator Instructions. 

3. System Development Aids 

The program written for the MCS-4 is metal-masked 
into the ROMs. Before doing this, the program needs to be 
tested to insure that it functions correctly in all situa- 
tions. To aid in testing programs written for the MCS-4, 
there are available three system development aids. These 
aids include an assembler and interpreter, a complete hardware 
prototyping system, and a resident software monitor. 

The MCS-4 Assembler and Interpreter is an ALGOL-W 
program which can be used to test and debug programs fOr 
the MCS-4 microcomputer. The system consists of an assem- 
bler which allows symbolic vrogramming of the MCS-4 ROM 
and an interpreter which simulates the actions of the MCS-4. 
The program provides extensive diagnostic Facalities. Lor 


monitoring the actions of the MCS-' program. 
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(Those instructions preceded by an asterisk (*) are 2 word instructions that occupy 2 successive locations in ROM] 
MACHINE INSTRUCTIONS (Logic 1 = Low Voltage = Negative Voltage; Logic 0 = High Voltage = Ground } 


OPR OPA 
0305 0, Dy D3 D, 0, Dy OESCRIPTION OF OPERATION 


JON 0001 jump to ROM address A2 A2 Az Ad, Ay Ay Ay Ay {within the same 
ROM that contains this JCN instruction) af condition Cy C2 C3 C4ll? 












Az Ao Az A2 is true, otherwise skip {go to the next Instruction In sequence}, 
0010 RRRO Fetch Immediate {direct} from ROM Data 02, 0; to index register pair 
Dz 07 03 0 0, 0,0, 0, location RRR, f2) 


Send register contro!. Send the address (contents of index register pair RRR) 
to ROM and RAM at X2 and X34 time in the Instruction Cycle. 

R RRO Fetch indirect from ROM, Send contents of index register pate location O 
Out as an address. Data fetched is placed Into register pair location RRR, 


JIN Jump indirect. Send contents of register pair RRR out as an address 
at Ay and Az time in the Instruction Cycle. 


An Ana AQAA 
*JUN o 100 33 633 Jui uNUNditiunal (Oo HUM address AZ AZ, A 
wsuiia ° 1- 
Az 42 Az Ag Ay Ay Ay Ay 
1 Oo 1 

















* IMS 0 AxzAzhg Ay Jump to subroutine ROM address Aa, Ap, Aj, save old address. (Up 1 level 
Az Az As Az Ay Ay Ay Ay In stack.} 
Mies i Sore | RRRR Increment contents of register RRAR, {3) 
ie ae RRRR Increment contents of register RRAR, Go to ROM address Az, Ay 
Aan te Re RRA {within the same ROM that contains this 1SZ instruction) if result #0, 
Da 2 Ar Taare Otherwise skip (go to the next instruction In sequence), 
| AOD R RRR Acd contents of register RRRA to accumulator with carry. 


SUB 


Subtract contents of register RRRR to accumulator with borrow, 

Load contents of register RRRR to accumulator. 

Exchange contents of index register RRRA and accumulator, 

Granch back (down 1 level in stack) and load data ODDO to accumulator, 


0©o090 Load data ODOD to accumulator. 
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Table IIL. MCS-4 Machine Instructions 
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Wirete the contents of the accumutator into the previously selected 
RAM output port. {Oulput Lines) 

Verite the contents of the accumulator nto the previously selected 
ROM output port. (1/0 Lines} 

Write the contents of the accumulator into the previously selected 
half byte of read/write program memory {for use with 4008/4009 only) 
Write the contents of the accumulator into the previously selected 
RAM status character 0 

Write the contents of the accumulator into the previously selected 
RA? status character 1, 

Write the contents of the accumulator into the previously selected 
RANI status character 2 

Viewte the contents of the accumulator into the previously selected 


RAN eestor charzevor 9 





Subtract the previously selected RAM main memory character from 
accumulator with borrow 

Read the previously selected RAM main memory character 

into the accumulator, 

Read the contents ol the previously selected ROM input port 

into the accumulator, (1/O tines) 

Add the previously selected RAM main memory character to 
accumulator with carry. 


Read the previously selected RAM status character 0 into accumulator, 
Read the previously selected RAM status character 1 into accumulator, 
Read the previously sctected RAM status character 2 into accumulator, 


Read the previously selected RAM status character 3 into accumulator, 


S-4 Input/Output Instructions 


Ny 








ACCUMULATOR GROUP INSTRUCTIONS 


1 1 1°43 010 1 


aa ere 


00 (0 (0 


RAL 


Oo 1 


DCL 


NOTES (Whe condition code Is assigned at follows 
Cy °1 
Cy *Q 


Invert jump condition 
Not Invert sump condition 





Col 
G,°! 


Clear both. (Accumulator and carry} 


increment accumulator, 
Complement carry. 


Complement accumulator. 





Rotate left. (Accumulator and carry) 

Rotate right. (Accumutator and catry) 
Transmit carry to accumutator and clear carry. 
Decrement accumulator 

Transfer carry subtract and clear carry. 

Set carry. 


Decimal adjust accumulator. 


Keyboard process Converts the contents of the accumulator froma 
one out of lour code to a binary code. . 


Designate command line, 


Jump If accumulator ig 2er0 
Jump If carry/iink is a 1 


C,4°* 1 Jump If test signal 1s 8 0 


(2}RAR is the address of 1 Of 8 index register palrs In the CPU. 


IJIR RAR is the address Of 1 of 16 Index registers in the CPU, 


{QE ach RAM chip has 4 registers, each with twenty 4 


-BDit characters subdivided Into 16 main memory characters and 4 status choractert. 


Chip number, RAM register and msin memory character are addreseed by an SRC Instruction, For the selected chip and register, however, 
status character locations are wlected by the instruction code (OPA}. 


Table V. 


MCS-~4 Accumulator Instructions 





There are two hardware prototyping systems avail- 
able. Each system uses the electrically programmable and 
erasable ROM (PROM). The 1701 or 1702 PROM is programmed on 
the prototyping system and simulates the exact action of 
a programmed ROM. The SIM4-02 prototyping board consists 
of a CPU and is capable of programming and testing up to 
sixteen 1701 or 1702 PROMs. The 4004 CPU also controls 
sixteen RAMs together with TTL simulations of eight ROM 
output ports and eight ROM input ports. The SIM4-01 is 
designed as a prototype for small systems. Tiers Hoard 
contains provisions for up to four le (OdlsmO r Olinda, Oe S - 

It also provides up to four RAM output ports, four ROM 
ourput ports and four ROM input ports. For small quantity 
systems where the cost of designing a metal-masked ROM can 
not be justified, the 1701 or 1702 PROM can be used instead. 

The recently announced INTELLEC-4 system makes avail- 
able to the programmer a resident software monitor, assem— 
bler, PROM programming, and text editing @apabilicy. 

4. Costs 

The cost of developing a system that uses the MCS—4 
Microcomputer can be divided into two areas: Hardware and 
software. The hardware costs have been shown to be small 
for the microcomputer. This is a list of the MCS-4 hardware 


costs: 
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1) (over 100) (over 1000) 


CPU $60 $30 each $15 each 
ROM $60 $15 each $ 5 each 
RAM $ 30 $15 each $ 5 each 


The software costsof programming the microcomputer 

are not really known. The assembler type language used to 
“<—program the MCS-4 requires the:programmer to keep track of 

the contents of all registers as the program steps through 
its set of instructions. Even one bit, such as the carry 
bit, can cause large errors to be produced it MOuweaccOuntbed 
for. All indexing and transferring of data between the 
processor and the assigned locations in memory must be 
feitten into the program. 

The task of writing a small program on the MCS-4 
is fairly straight forward, however, as the Compu vavlonat 
complexity grows, the task of writing the program increases 
mapLldaly. 

The sections following this discuss the problems 
encountered in programming the NCS-iwto menales the computa 
tional tasks of a large navigation system. Aids in software 


development were investigated and are also covered. 
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III. NAVIGATION SYSTEM 


Air Navigation is the process of directing the movement 
of an aircraft from an initial point to a desired final 
moing, A coordinate reference frame must be established 
and the initial point and the final point must be located 
in the reference coordinates. A Navigation system must 
provide timely coordinate measurement and computation of 
fecirea aircraft position. 

The prime navigation problem years ago was to reach the 
desired destination. Today many aircraft wish to use the 
same airspace at the same time. This problem has required 
Peperavion svandarus to ve put invo efrecc. Pi gon. 
Msciemed a block Of space, measured by laveral, loneiGudinmailk, 
and vertical dimensions, which moves at the speed indicated 
Mmuneipilot's flight plan. It is the pilot's responsibility 
£0 remain within this block of space. As more aircraft 
wish to use the same airspace, the separation standards will 
have to be reduced. The desired navigation system is one 
eat eBives a continuous, real time indication of where the 
aircraft is located. Navigation systems must be required 
to have greater accuracy, greater automation, and simplicity 
in operation and display. 


Dieenecanwor a navigntaon System to be reliable is a 


necessity. Equipment failure must not endanger the ay rerar us 


ig 





The reliability of a navigation system can not be measured 
solely by the fault-free operation of one piece of equipment. 
Reliability is a function of the operation of the total navi- 
Beaeton system. In case of partial equipment failure, the 
navigation system must automatically switch to an alternate 
source of information. Therefore, decision making el reulLery 
must also be a part of the navigation system to insure @2 
continuous flow of accurate navigation information. 

The following is a review of current navigation systems 
available. From this information, the possibility of inte- 
grating the MCS-4 microcomputer into a total navigation 


system will be investigated. 


A. SELF-CONTAINED NAVIGATION SYSTEM 

Only automated navigation systems can Savicty coewkeduuce— 
ment for continuous and accurate navigation information. 
The need is for a pictorial type of display which will give 
the pilot an accurate and immediate indication of the air- 
craft's present position. Economic constraints have en- 
couraged that this navigation system be standard throughout 
the military and have a close commonality with Civ teaovi alo 
needs. Many of the military and civil aircraft in use today 
have space provided for this type of system. 

The navifation systems in greatest use today are PaG.o 
Mavigation syst@éms. In order to keep picmcoct ana complexivy 
of airborne radio navigation equipment down, a comprehensi.ve 


system of large and costly radio stations must be kept on 


Hg 





the ground. Some of the major systems in use today are 
listed in Table VI. When such a system is created, it is 
almost impossible to abolish it. The high iabaal cost: OL 
setting up a new system and the large number of aircraft 
equipped for and dependent on the system in use, make it 
hard for new systems based on new technology to become 
established. 

Radio navigation systems have many limitations. Some 
of these limitations are due to the short distance On eea-ouncd 
wave propagation, sky wave contamination, atmospheric noise, 
Bulltaoath effects, and site error which can result in ambi- 
puious position fixes. Errors in radio navigation signals 
can not be predicted because they are a fanetuom Of dally 
ana seasonal changes in environment, temperature, ionosphere 
ecation, and local weather. Even the placement and type 
of radio used in the system is affected by political factors. 
These limitations, together with frequency interference 
between radio stations and the susceptability of radio 
navigation systems to jamming, all support the desirability 
of aircraft to have a self-contained navigation system. 

Military aircraft navigation systems require worid-wide 
flexibility with the ability to navigate Vane MOUG ce uLanee 
on ground-based aids or the use Of equipment susceptible 
fomamming. Civil aircraft that fly over oceanic or desert 
routes require a self-contained navigation system as ies 


primary source of position fixing. The development of 
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hovercraft that must cross large bodies of water without 
land based navigation aids and the need for cruise missiles 
that can independently maneuver over hostile territor all 


point to the need for self-contained navigation systems. 


By ALR DATA SYSTEM 

An air-data system consists of aerodynamic and thermo- 
dynamic sensors inputting to a central air-data computer. 

The sensors measure the characteristics of the air surrounding 
the vehicle and input them to the computer. The computer 
calculates flight parameters such as true airspeed, free- 
stream outside-air temperature, and Mach number. The sensors 
required in this system are angle-of-attack vane, static 
Mesure source, pitot tube, and total temperature probe. 

In order to use tne air-data computations for dead- 
meckOning, the attitude and heading of the vehicle must be 
supplied from some external source onboard. This can be in 
the form of a simple directional and vertical gyroscope or 
a stable-platiorm configuration such as an Inertial Navigation 
System. TALseiIntormatten tOgetner With information on the 
velocity and direction of the wind allows the navigator to 
determine the position of the vehicle by extrapolating from 


meoreviously known fixed position. 


fm DOPPLER NAVIGATION SYSTEM 
A Doppler Navigation system radiates a pattern of beams 


to the surface of the earth and receives the reflection of 





t 


this energy back. The difference between the frequency of 


the signals transmitted and the frequency of the signals 


received is called the Doppler Effect and can be used to 


compute the vehicle's velocity along each beam. The 


frequency shift in each of the beams is detected and used 


by the system computer to calculate the distance traveled 


along and across the vehicle's true heading. 


The advantages of the Doppler Navigation System, as 


listed by Kayton and Fried [Ref. 21], are as follows: 


ae. 


Oi 


) 


) 


The 


lie 


oe 
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It provides continuous velocity and position with 
respect to the ground. 


It is completely self contained. 


“Its average-velocity information is extremely 


accurate. 
Doppler-Radar information is obtainable anywhere 
on earth, including over oceans. 


It 1s an all-weather system. 


Doppler Radars are amenable to high-reliability 
all-solid-state design because of their low 
radiated power. 


The Doppler Navigator Goes not require preflight 
alignment or warmup. 


It radiates at microwave frequencies. 


i 


disédvantages of the Doppler Navigation System are: 


ft is dependent on an external direction sensor 
fOr soe imueh Ln ermal ron. 


Ene  WOSition Jnioumavion derived from a Doppler 
Navigator degrades as the distance traveled 
increases. 


The short-term or instantaneous velocity information 
is not as accurate as the average velocity. 





D. INERTIAL NAVIGATION SYSTEM 

Airborne inertial navigation systems, hereafter re- 
ferred to as INS, operate on the principle that every time 
the vehicle changes speed or direction it is said to experi- 
ence an acceleration. By measuring this acceleration, the 
velocity and distance traveled by the aircraft can be found 
by integrations. 

The INS is a gyro-stabilized platform with accelerometers 
mounted so that accelerations of the vehicle are measured 
in the north/south, east/west and vertical directions. The 
total system is mounted on gimbals to allow the vehicle to 
rotate without disturbing the attitude of the stable platform. 

The INS must be kep tangent to the earth's surface at 
all times or the accelerometers will experience an acceiera-~ 
Gon error due to gravity. This problem is solved by con- 
Structing a theoretical pendulum in the system with its 
POvmcieeauniemeCnlci Ol Une cartun and tts Oener end on the 
Surface of the earth. The aligning of the INS before flight 
wall bring this pendulum to rest making it impossible to 
be set in motion by any force at the surface of the earth. 
Any small error in the alignment causes a sinusoidal error, 
called the Schuler Pendulum effect to be introduced in all 
imertials. This effect is periodic with a period of 84 
iMenutes. 

The Inertial Navigation System has the following advan- 


tages as listed by Kayton and Fried [Ref. 21]: 





a1;) its indications of position and velocity are 
instantaneous and continuous. 


2.) It is completely self-contained. 
3.) It is nonradiating and nonjammable. 


4.) Navigation information is obtainable anywhere on 
Dre -eanun. 


5.) It is an all-weather system. 


«6.) -Navigation-information is -independent of the 
vehicle's maneuvers. 


7.) It directly provides outputs of position, ground- 
speed, and vertical position. 


SO eitelts the Mostuaccurate source of roll. pitch, 
and attitude of the vehicle. 


The Inertial has the following disadvantages: 


i 
Cy 
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Pete pescitioOm ano vetourly silicotnacicin Gerraces 


with time. 


2.) The equipment is expensive and relatively difficult 
to maintain and service. 


Bo) The INS muss De anitially aligned. 


E. POSITION FIXING SYSTEMS 

Position fixing is the determination of the vehicle's 
Pestti1on without reference to any former position. short 
range radio systems such as ADF, TACAN, DECCA, VOR, and 
LORAN A are not considered applicable for long range 
navigation. 

1. Long Range Radio Navigation Systems 

The two primary long range Radio Navigation Systems 


in use today are LORAN C and OMEGA. LORAN C and OMEGA are 





both hyperbolic fixing systems. The vehicle's fix position 
is obtained by measuring the relative distance from two or 
more stations of known ae 
a. The LORAN C Network 

The LORAN C network is comprised of master ground 
station and at least two associated slave ground stations. 
Pmeeea transmissions radiated=from the master station are 
received at the slave stations; each slave station then 
transmits similar groups of pulses, synchronized accurately 
with the signals received from the master station. Receivers 
on-board the vehicle measure the time differences between 
the master and slave's transmissions which allows an automatic 
ime COMDULaba ON . 

The following are advantages and capabilities 
of LORAN C: 


1.) Currently operational with minaturized receivers 
developed. 


Pe) ee hence tom! s a0 Mm. 


3.) Continuous fixing provided. 


4.) Passive at the receivers. 

5.) All weather capability 

Ge eccimece oOSteilOon 11% to ;cO0Tt. 

7.) Total cost of on-board receiving equipment 


approximately $10,000. 


The disadvantages of JLORAN C are: 
1.) Position fix in hyperbolic coordinates. 


2.) Dependent on fixed pfround stations. 
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3.) Susceptible to atmospheric noise and sky-wave 
contamination. 


4.) Susceptible to jamming. 


5.) Total world coverage ‘not available. 


b. The OMEGA System 

OMEGA is an earth reference hyperbolic phase- 
Matching navigation system operating at 10.2, 11.33, and 
13.6 kHz. The system is designed to provide world-wide 
coverage with accuracies better than one nautical mile. 
The position of the vehicle is determined by measuring the 
relative phase of the Signals transmitted from two or more 
Stations. 

OMEGA has the following advantages: 

Range 8000 nm. 
Coverage is global with just eight stations. 


) 
) 
So meconmtinucus taxing avarlaple, 
) Passive receiver on the vehicle. 


5.) All weather capability. 


The disadvantages of OMEGA are: 


te eerOstticon fix is in hyperbolic coordinates. 


NO 
we 


Dependent on fixed ground stations. 
3.) Susceptible to Ionospheric disturbances. 
4 


-) Susceptible to jamming. 


e. The Satellite Navigation System 
oatellite havigation is a method of fixing the 


vehicle's position from data obtained from an artificial 





earth satellite. The position of an artificial satellite 
can be predicted by the orbit it is traveling. A relative 
fix can be made from the satellite by measuring the elevation 
angle, the azimuth angie, range to the satellite, or the 
rate of change of altitude, azimuth or range. If a single 
variable is measured, Successive measurements are required. 
There are three basic types of satellite navigation: 
Doppler Systems, which measure rate of change of slant range; 
angle measuring systems, which measure elevation angle; and 
manging systems, which measure the slant range. 
Satellite Navigation Systems have the following 
advantages: 
1.) All-weather position fixing. 
2.) Communication service to the vehicle while in flight. 
3.) Worldwide service. 
mT capable ot accuracies up to 200 ft, 


5.) Can provide vehicle with velocity and heading 
information as well as position information. 


6.) Passive at the receiver. 


The disadvantages of a Satellite Navigation System 
are as follows: 


1.) Special receiving equipment required including 
antennas with vertical patterns. 


2.) iiescacellige orbit must be predicted and 
continuously up-dated. 


3.) The satellite's signal may be susceptible to 
jamming. 


4.) The satellite signal is susceptible to Ionospheric 
retraction. 
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5.) Navigation errors are introduced by drifts in the 
frequencies of the satellite transmitters. 
3. The Terrain Mapping System 
A Terrain Mapping System operates on the principle 
that the geographic location of any place on the land surface 
of the earth is uniquely defined by the vertical contours 
or voporrapny of “the*surrounding ~area. “Phe Terrain Mapping 
system measured the vertical contour of the terrain along 
its flight path, using a radar altimeter to measSure clearance 
ieeve tic cteorrain end a reference altitude of the vehicle. 
By subtracting instantaneous radar measured altitude from 
the reference altitude, the Terrain Mapping System determines 
the terrain contour. The system then searches its computer 
memory to find a stored terrain contour, whose coordinates 
are known, which closely matches the measured one. This 
serves to fix the vehicle's position. 
The advantages of the Terrain Mapping System are 

as follows: 

1.) It is an all-weather system. 
.) It is completely self-contained. 


Ties snereibLyeareslStant tO Jamming. 


ag a Se, 
~~ 


Pesltetsvoccurate wip. to. 50) ft. 


Accuracy demands imposed on the radar altimeter are 
not excessive. 
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6.) It is unaffected by man-made changes in topography. 


The Terrain Mapping System has the following 


disadvantages: 


1.) The terrain contour to be flown over must be 
previously mapped. 


2.) To increase the accuracy of the system, there must 
be an increase in the computer memory Gapacity. 


3.) <A dead-reckoning system must also be carried so the 
Terrain Mapping System knows where to search for a 
seeontour*™match. 

4.) It can not be used to obtain position fixes over 
the ocean or large bodies of water. 

me THE INTEGRATED SYSTEM 

The basic parameters of the navigation systems currently 

available have been reviewed. In order to select the proper 


long range navigation system for military use, a minimum 


set of system requirements was established: 


Range 3000 nm. 

Accuracy Soin. 

Velocity O- = 2000 kts: 

Coverage Worldwide 

Weather Restriction All-weather capability 

Pees penal y Capable of self-testing and 
switching mode when needed 

Equipment Non-jammable, completely 


self-contained 


A review of the previously mentioned navigation systems 
indicates that no single system available is capable or 
flexible enough to meet the system requirements. The Doppler 
Navigation System is self-contained and provides very accu- 
rate average velocity measurements, however, the system can 
not sense short term velocity fluctuations and must depend 


on an external heading source. The Inertial Naviration 
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System is self-contained and provides excellent short- 
term position, velocity, and heading information. The 
Inertial Navigation System, however, has a long-term error 
build-up due to the inherent Schuler Cycle in the system. 

In order to meet the system requirements, an integrated- 
navigation-system was chosen. An integrated-navigation- 
system feeds the outputs of several navigation systems into 
a central computer which then provides a single more accu- 
rate output. The central computer provides decision making 
ability to test the accuracy of the individual navigation 
systems. By inputting fixed positions to the central com- 
puter at regular intervals, the computer can determine 
system drift varameters and provide the proper system bias 
Pelocity. 

The navigation system developed in this paper is an 
Inertial/Doppler system integrated by an MCS-4 microcomputer. 
The microcomputer combines the short-term accuracy of the 
Inertial with the long-term accuracy of the Doppler to 
eptain the most probable position of the vehicle. The 
micro-computer is programmed to provide decision making 
flexibility to ensure the outputs remain accurate during 
Partial system failure. The wind influencing the vehicle 
is continually computed and updated in the MCS-4 memory so 
that the computer can automatically switch to an air-data 


mode Of operation in case of Inertial and Doppler failure. 
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The navigation equations used to compute the vehicle's 
position from the inputted data are covered in the next 


section. 
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IV. NAVIGATION EQUATIONS 


The navigation equations necessary to calculate the 
position of the vehicle relative to the Bait are presented 
in this section. The mathematical model presented takes 
the dead-reckoning outputs of the Inertial, Doppler, and 
Air-Mass Systems and extrapolates the present position of 
the vehicle from the last known position. Direct. position 
data can be inputted into the algorithm and is used to update 
the position of the vehicle. The system drift error is 
computed by comparing the known position with the dead- 
reckoned computed position. 

meee cCh re ime meuoGlcoLoOns 1m Gevcll. @ Sra0hical oiec- 
ture of each equation's function is included. The graphical 
picture represents each equation as a black box with speci- 
fed inputs which produce the deSired outputs. A complete 
meseription of graph theory as applied to this paper is 


given in Section V-A. 


A. THE EARTH MODEL AND THE COORDINATE REFERENCE FRAME 
inewiarst. Step Imnvdeveloping the Navigation equations 

is to prepare a mathematical model of the earth. The earth 

can be approximated as a sphere of radius R, the nominal 

mwa vOrlial radius. However, this approximation would produce 

mm) waacceptable error as great as 3 miles for every 1,000 

Mieeureveiread.,. Ihis 1S due to the flattening of the earth 


at the poles and the bulging of the carth at the equator. 
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A more exact model of the earth is the reference 
spheroid where the earth is assumed to be an ellipsoid 
of revolution with the semi-major axis R, the nominal 
equatorial radius, and the semi-minor axis of radius P, 
ene polar radius. Accuracy to within 30 ft. for every 
mOOO nm. traveled is possible with this model. MThe 
determination-of angels, however, in this system is very 
set ficult. 

The model chosen in this system is a combination of the 
lwo previously stated models. All angles are determined 
from the spherical model thus allowing the convenient lati- 
tude/longitude coordinate reference frame to be used to 
map each point on the real earth to the earth model. The 
reference spheroid is then used to determine arc lengths 
that more closely represent the actual shape of the earth. 
The arc lengths are applied to the determined latitude and 
fongitude in the form of conversion constants. This method 
of computation is very simple and has a maximum error of 


mo ft. tor every 1000 nm. traveled. 


B. DERIVATION OF CONVERSION CONSTANTS 

tiem mueriauuona i nauulcolemule equals 6076.1033 ft. 
One minute of arc length measured at the equator equals 
6087.08 ft. The longitude conversion constant, KLO, is used 
to relate nautical measurement to the actual geographical 


measurement: 
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KLo = Leet/degree 
meaty Nautica le mile 


KLO = 6087.08 feet per minute/60 minutes per degree 
076.1033 feet/l nautical mile 


60.1084 nautical miles/degree 


I 


KLO 


The average length of a minute of arc length on a 
meridian is 6076.82 ft. To determine the latitude conver- 
Si0n constant, KLA, which 1s also used to relate nautical 
measurement to the actual geographical measurement, the 


following equation is used: 


. feet/degree 
KLA fest Nautical mile 


6076.82 feet per minute/60 minutes per degree 


BLA = COmoOnemmce ty ll maniuea le mame 
KLA = 60.0071 nautical miles/degree 


iWewdtstance traveled in nautical miles, east-west and 
north-south, is divided by the conversion constants resulting 
in distance traveled measurements in degrees, longitude and 


feacitude. 


C. WIND CALCULATIONS 
The wind effecting the movlon of the vehicle can be 


obtained by comparing the north-south and east-west vector 
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components of the true airspeed, as computed by the air- 
mass system, with the vector components of ground speed, as 
computed by the inertial system. The calculated wind com- 
ponents mist be continually updated so the current wind 
components stored in the computers memory are the most cur- 
rent. In case of failure of the Inertial and Doppler Sys- 
tems, the computer will be capable of computing the dead- 
reckoning position from the outputs of the air-mass system 
and the last stored wind components. 

The wind velocity components are computed from the 


following equations: 


TAS VAX WIND 





Components VAY VELOCITY Vik == 
INERTIAL Meee COMPUTATION VWwik -———- 





I 


VWXR VGXI -— VAX 


VWYR = VGYI -— VAY 


The wind velocity components computed are subject to 
musts, Maneuvers, and other short term disturbances. In 
eraer to decrease the effects of random fluctuations in the 
wind velocity, the new computed wind velocity is averaged 


Mith the previously computed wind velocity in a smoothing 
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routine. The smooth wind velocity is derived from the 


following equations: 








SMOOTH 
WIND 


—— VWYR 
—WiX, 5 VELOCITY 
computation | VW¥s-—> 
ey 
i-1 

aa 

VWX = G VWX + (1-G) VWXR 

VWY = G VWY + (1-G) VWYR 


D. DEAD-RECKONING CALCULATION 

Dead-reckoning is a means of navigation in the absence 
of position fixes. The vehicle's position can be estimated 
by measurements of the groundspeed components. The ground- 
speed components are integrated over a given time interval 
to give the distance traveled in that time interval. The 
estimated position of the vehicle is the Summation of the 
agistance increments traveled from the last known fixed 
Bosition. 

1. Air Mass Mode 

The air mass mode of operation is a back up for the 

two primary dead-reckoning navigation systems, the Inertial 
end Doppier. The components of airspeed, true heading, and 
stored wind information are combined to output distance 


increments in the north-south and east-west directions. 
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The true air-speed inputs are influenced by random 
fluctuations in the Pitot-Static System. In order not to 
influence the navigation cycle by these fluctuations, the 
inputted true air-speed is averaged with the previously 
stored true air-speed. The smooth true air-speed is derived 


from the following equation: 






Aen 
—— TAS. 
i= 


SMOOTH TRUE 
AIR-SPEED 
CALCULATION 







TAS 


ASS a eo 


TAS = 2 


The north-south and east-west components of air-speed 
are calculated using the previously calculated true-air-speed 
and the true heading inputted from the INS. In case of INS 
failure, magnetic heading may be inputted from the flux 
valve together with the magnetic variations applicable set 
in by hand. The air-speed components are calculated as 


follows: 









COMPONENT 
AIR-SPEED 


CALCULATION 





Ve eer 


i —— 





VAX 


TAS SIN (TH) 


VAY = TAS COS (TH) 
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The air-speed components represent the velocity and 
direction of the vehicle through the air-mass. The ground 
track of the vehicle can be found by adding vectorally mae 
motion of the air-mass, wind velocity, to the air-speed as 
shown in Figure 3. The equations for the ground velocity 


tn the air-mass mode are as follows: 
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VGXW = VWY + VAX 


VWY + VAY 


i 


VGYW 


Lastly, the ground distance traveled by the vehicle 
during one navigation cycle of the computer can be determined 
by integrating the ground velocity over the time interval of 


eae navigation cycle. 
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Figure 3. NAVIGATION VECTOR DIAGRAM 
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ee lcepter—Mode 

The Doppler System obtains very accurate average 
velocity measurements and converts them to distance measure- 
ments along and across the vehicle's true heading. These 
measurements along with true heading, from an external 
meurce, can be used to calculate the distance traveled by 
the vehicle as shown in Figure 3. 

The dead reckoning position increments computed from 


the doppler are derived from the following equations: 


Doppler SDD 
SHD 


Inertial Te 













Doppler 
Distance 
Computation 


DXD SP PmocueOoe( TH )) + chp x OPH) 


DYD SD emo Nl) te Des N= CEH) 
3. Inertial Mode 

The Inertial Navigation System, (INS), provides 
excellent short-term velocity and heading information. The 
outputs of the INS are the velocity components of the vehicle 
north-south and east-west. The ground track of the vehicle 
can be determined as shown in Figure 3. 

The dead reckoningg position increments computed 


mom the inertial are derived Srom the following equations: 


TO 
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Distance 
Computation 
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4, Integrated System Mode 
To take advantage of the extremely accurate average- 
distance measurement of the Doppler and the-precise short- 
term distance measurement of the Inertial, the distances 
measurement by each is averaged in a smoothing routine to 
give the most accurate dead-reckoning distance measurements 
possible. The smoothing factor, G, used in the smoothing 


equation is set depending on the accuracy of each system. 
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E. POSITION CALCULATION 
iPr ocOg@ophie Hix FOsitivon 

Any time the vehicle's position is fixed, the updated 
fixed position will be loaded into the computer's memory 
along with the time the fix was taken. The dead-reckoning 
position coordinates at the time of the fix will be up-dated 
to the new position coordinates and the computer will continue 
to dead-reckon the vehicle's position from the last inputted 


fe xX. 










Update 
Nav 


LOF Position 


—— LAF 


LON ————> 
LAT ———> 






LON = LOF 


Cal = LAP 


Ewe Up. Wave ceographic Fost tion. 

Thies anevea poste Cn set wine -venwcle 15° loaded into 
fie computer memory before launch. The geographic position 
is up-dated each time a fix is taken. The geographic posi- 
tion between fixes is updated by dead-reckoning distance 
computations from the Doppler and Inertial. The zero sub- 
ferupt indicates the vehicle's previous geofraphic position. 
The computed system drift increments are also inputted to 
account for system drift crrors in the dead-reckoning 


computation. 
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DX + SDX 


LON = LON, + ero-x 0S (TIAT) 


Sa wovyStem Drift Error 
Each time a new position fix 1s taken, the fixed 
position is compared with the calculated dead-reckoned 


Mesition. From this information, the amount of drift in 


the dead-recKnxoning computations 


=O 
LAT Dray SDYC -————-> 
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tea ovem olas, Velocity 


The system drift error computed from the previous 
mex 15 used to update the system bias velocity. The system 


Mees Velocity iS used to cancel the system drift error. 
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5S. System Bias Distance 
The system bias distance is computed by integrating 
the system bias velocity over the time of one nav cycle. 
The system bias distance is applied to the up-date geographic 


position to cancel the effects of the system drift. 
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The navigation program combines these equations as 


described in the Navigation Functional Flow Chart, Figure 5. 
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V. PROGRAMMING THE MICROCOMPUTER FOR NAVIGATION 


‘The Navigation equations discussed in Section IV des- 
Ceribed how to process the outputs of the Inertial, Doppler, 
and Air-Mass Systems in order to calculate the current 
position of the vehicle. In order to make a working system, 
the Navigation equations had to be programmed into the 
software of the MCS-4 microcomputer. 

The program was designed to incorporate the following 
principles and requirements: 


1.) The total program had to be written in the assembly 
language of the MCS-4 microcomputer. 


nat 
4a’ VV 


2.) The total run time of one navigation cycle could 
exceed 200 milliseconds. 


3.) The design of the program should be modular to 
facilitate the addition or deletion of new code. 


4.) The program should consist of an executive routine 
that calls upon various subroutines as indicated 
on a process graph of the navigation program. 

be) her desien, allows ior SubSsticution or addition 
of new functions, routines, and/or programs in 
a Straight forward manner. 

6.) For each built-in function, a single subroutine 
is called to perform the calculation. 

The instruction set of the MCS-4 microcomputer, as des- 
erabed in Section J1I.C., contains only addition and subtrac- 
tion of two four digit binary numbers as the fundamental 
arithmetic operation. In order to make the detailed pro- 


gramming easier, a so-called "process graph" method was 


maea to cevelop the prograin. 





Mm GRAPH THEORY 

Programming the navigation equations would be an easy 
task in a higher level language such as FORTRAN. To program 
the same equations in an assembly language requires a detailed 
knowledge of the current contents of all index registers 
and memory locations. The status of even one Single bit 
suchas "the *carry~bit can not-be left unaccounted. The 
requirements of keeping minute details in mind led to the 
development of a graphical means of representing the 
microcomputer program. 

Graph theory provides a Simple and powerful tool for 
constructing mathematical models of discrete arrangements 
of objects. The process graph consists of vertices which 
are pairwise connected by a directed line. 

The best means to illustrate the concept of graph theory 
applied to microprogramming is through an example. Consider 
the problem of adding two numbers together which reside at 
Biven locations in memory and whose sum is to be located 
in a specified location in memory depending on its sign. 
Mats process can be described by the following flowchart: 


No 
Lo tL Leo. ——— > 


Whas notation does not give much information about how the 


Processor 1S to™waccomplish this process, how much memory is 





required by the program, the input and output variables, 
and how much time is needed to carry out the process. 
The same example written aS a process graph could be 


as follows: 


Re et 
Z, 
X 
ip 
y ae _ 
Z, 


The process graph indicates that the variables X and Y are 
inputs to the function add and Z is the output of that func- 
tion. A test is made of the output, Z, to decide where Z 
mea we locacea. The amcounc ci time and memory involved in 
this computation can be easily established by referencing 
the known memory space and time required for each operation 
mavOolved. The depth to which the programmer draws the re- 
Suired graph for his computation is a function of the capa- 
bility of the operations he has available to him together 
with his own preferences. The same example discussed could 


mee fraphed in more depth, as follows: 


RAM. 
ACC ee 
Acc RAM 30 





a 





Another programmer reading this graph only needs to know 


the following in order to write the required program: 


ACC 





+ 
ACC 


i: RAM 
ae 
ACC 


Load the variable X 
into the accumulator 


Add the variable Y to 


the contents of the 
accumulator 


Mest the Siegen of the 
contents of the 
accumulator 


it the Sven C1 the 
Aecumubaror 1S DGsili vee 
store the contents of 
the accumulator into 
RAM location (2,0) 


if ene Sian. of the 
accumulator is negative, 
store the contents of 
the accumulator into 
RAM location (3,0). 


The detailed process graphs used in developing the navi- 


gation program were a great aid in the design, analysis and 


documentation of the program. 


The following is a list of 





the problems that could be analysed directly from the 
process graph: 
1.) What portions of the problem are sequential. 


2.) What portions of the problem may be processed in 
parallel by independently Operavlnie proce sSOrS: 


3.) Where and how intercommunication may take place 
when a multiplicity of processors is used. 


4.) What capacity must the processor possess before it 
can be successfully used to solve the problem. 


5.) What is the maximum allowable time for the operation 
of each subprogram. 


6.) What is the memory space requirement for each 
Suu Ou elle. 
B. DEVELOPMENT OF NAVIGATION PROCESS GRAPHS 

The MCS-4 microcomputer was first thought of as a black 
box taking inputs from the Inertial, Doppler, Air-—-Mass, and 
Position Fixing Systems and outputting the vehicles Clee 
position in latitude and longitude. The Pe a) eC Ge lae ae OileG moles 
tation of this system is showm in Figure 4. 

Figure 4 represents the total navigation system. The 
Inertial, Doppler, Air-Mass, and Position Fixing Systems 
mrovide input to the system as described im soeICu LOM LV. 

The microcomputer program was required to complete the sys- 
tem. The programming was accomplished by continually break- 
ing down the functional operations Te omeicadle and smaller 
parts until the operations were simple enough to be directly 


written in the MCS-4 machine language. The functiona] process 
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graph, Figure 5, represents how the microcomputer program 
was initially broken down into the fundamental navigation 
equations es described in Section V. 

The functional process graph was analyzed to define the 
flow of the input variables and to determine the feasibility 
of a multi-processor system in order to shorten the required 
computational time. It was noted that the program could be 
broken into two parts, the calculation of the distance in- 
erements, and the calculation of the latitude and longitude 
from the distance increments. The analysis demonstrated 
that the computation time could be nearly halved by having 
two microcomputers working simultaneously to produce the 
Gesired outputs. The system was designed to have one micro- 
computer receive the given inputs and compute the distance 
increments traveled, while at the same time, the second 
microcomputer computes the latitude and longitude from the 
previously calculated distance increments. 

The functional process graph was then broken down into 
the operational process graphs for each microcomputer, 
Figure 6, and Figure 7. The operational process graphs 
represented the desired program for each microcomputer. 

The program described by the operational process graph, 
Figure 6, was written to investigate the time and effort 
required to develop the required software. 

The program wes written in a modular form as required. 


Each operation represented in the operational process graph 
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was written as a separate subroutine. An executive routine 
was then developed to call each Subroutine in the proper 
order and account for each variable as required by the 


operational process graph. 


C. PROGRAM ANALYSIS 

A program analysis was developed to define those pro- 
blem areas that had to be solved before programming the 
microcomputer. The areas that were investigated were com- 
putational speed, memory space available, and accuracy 
required. 

The operational process graph, Figure 6, represented 
the tasks to be accomplished. The type of operations and 


Mumder Of Operations that were required are listed as follows: 


Operations Times Called 
Multiply 12 
Cosine i 
Sine af 
Addition 2 
Subtraction 3 
Division bys Guo 7 

Total 36 


The operational process graph was used to determine the 
poeeed limitations on each operation programmed. The total 
navigation cycle was limited to 200 milliseconcs. Since the 
number of operations required to be performed differed depen- 


ding on what sensors were operational, the critical path of 





the operational process graph was determined. The critical 
path occurred when the Inertial and Doppler were both opera- 
tional, Figure 8. The operations involved in the critical 
path were ten multiplies, Cosine and Sine calculations, ten 
Bacay ltOns, three Subtractions, and Seven divisions by two: 
The total computational time for the critical path was 
ebtmibed t0-200emildiseconds. 

The critical operations which had to be developed were 
fae Multiply, Cosine, and Sine subroutines. Investigating 
previous work on the MCS-4 indicated that previously pro- 
grammed multiply, Cosine, and Sine routines were requiring 
50 msec., 650 msec., and 750 msec. respectively. In order 
to program the microcomputer for navigation, the development 


of the following routines were required with the following 


POnstvraints: 
Poa oy less than 200 msec. 
Cosine less than 200 msec. 
Sine less than 200 mSec. 


The time used by each of these subroutines determined the 
number of processors required to meet the 200 millisecond 
mime constraint. 

the amount Of memory avallable to program the navigation 
routine was a function of the 4004 CPU and the number of 
microprocessors used. One 4004 CPU can directly drive six- 


teen ROMs and sixteen RAMs. The number of instructions in 
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the navigation routine was limited by the space available 
in the ROMs. Since each ROM could hold 256 instructions, 
the program was limited to 4,096 instructions per micro- 
processor. 

The time required to sequentially execute every instruc-— 
tion in the sixteen ROMs would be 44 mlliseconds. Since the 
program was required to be written with an executive routine 
and a set of subroutines that would repeat the same set of 
instructions several times, the limiting time constraint 
would be reached before using up the available ROM space 
in one microprocessor. It was determined from this analysis 
that ROM space would not be a limiting factor in writing 
imme navigation program. 

The amount of memory space available to store the values 
of each variable was determined by the space available in 
the RAMs. Since one 4004 CPU could drive 5120 bits of RAM, 
the navigation program was limited to 320 variables of LG 
bits for each microprocessor. The number of variables re- 
Quired was Getermined from the operational process graph, 
Figure 6, where each line connecting a pair of vertices 
represents one variable. There were puedes) Indacatine 
that a maximum of 54 variables were required plus those 
wmariables used ve any single operation. By overlaying 
weriabpbes in the same RAM memory space, the memory space 
requirement was reduced. It was determined from Fhis analysis 


that RAM space available would not be a limiting factor. 
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The required accuracy of the navigation program was a 
function of the accuracy of the input variables. The accu- 
racies of Inertial, Doppler, and Air-Mass systems used on- 


board the P3C aircraft were used in this analysis as repre- 


senting the state-of-the-art systems in naval aircraft today. 


The specifications for these systems are as follows: 


system Designation Accuracy 

Inertial ASN-84 Teo enous 
True Heading ASN~84 + 9 ARC-MIN 
Doppler APN=187 he) chiens 
True Air Speed Pitot-Static fie. 0 nor s 


The accuracy of the navigation program was a function 
Se une aeeurecy Cf input data es well ab the bit size 
assigned to each variable. The limited accuracy of the 
input data permitted each variable to be no greater than 
16 bits. This allowed each variable to be represented by 
four hexadecimal-digits with the first bit assigned as the 
Sign bit. The hexadecimal point for speed measurements was 
fixed so that there is one hexadecimal digit to the right 
of the decimal point. This allowed the accuracy of the 
Peecad inpues 7oO be within + 2.0625 knots. The range on the 
inputs due to a 16 bit variable limitation was + 2047.99 
knots. The accuracy requirement was not considered a major 
fimication in the program anelysis. 

The program analysis pointed out the major areas that 


had to be demonstrated feasible before programming, of the 
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microcomputer for navigation was initiated. The program 
analysis was used to set the design goals at each step and 
the process graph proved to be a major tool in performing 


the program analysis. 


D. SUBROUTINES 
A large number of the mathematical operations required 
in the navigation program are repeated many times. In order 
pO decrease the total programming effort and also decrease 
the memory-capacity requirements, many of the operations 
required were written as subroutines. These routines con- 
sist of a series of instructions dedicated to a specific 
task. The subroutines developed for the navigation program 
PSamaimiGeaqsimeo uWwO ErOuce. tmcse anvelving complex methe= 
matical operations and those involving more common func- 
tional operations. Each subroutine was written in a general 
form to permit its recurring usSe by the executive program. 
Pee jen titel Cab On ROutine 

The major limitations of the MCS-4 microcomputer to 
be overcome were the limited Ae rie set and slow speed 
piaeccweonlatton. |lhe most powertvl aravthmetic instruction 
available was the single addition two Reece The mul- 
tiplication routine was written to increase the capability 
of the MCS-4 microcomputer in order to satisfy the require- 
ments Of the navigation program. 

A multiplication routine had been written involving 


multiplication by a seri@s of additions. The program required 
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only fifty instructions, however, the computational time 
was 40 milliseconds. Since the navigation program required 
a minimum of ten multiplications, this method was unsatisfactory. 

One of the advantages the MCS-4 microcomputer has 
is its inexpensive memory. It was decided to investigate 
a different way of programming the microcomputer that would 
«take advantage of available memory. It was discovered that 
memory space could be traded for speed by a table look-up. 
The multiplication subroutine was written using a table 
look-up scheme. 

a. Capability of Subroutine 

The multiplication subroutine was designed to 

take a four hex-digit number, X, and multiply it by a four 
hex-digit number, Y, resulting in a four hex-digit chopped 
number. The input X is loaded in Index Registers R8 + RB 
with the least significant digit, Xo, loaded in R8. The 
Input Y is loaded in Index Registers RC + RF with the least 
Ploiiit@cany adipic, 1.o, loaded in, RC. The product will appear 
imeimndex Registers RC > RF wath the least significant digit 
Peecareqd an RC. 

b. Description of the Routine 

The process was done by a table look-up scheme 

Simitar to tonghand multiplication. An’example of this 


method follows: 
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nr X x es) 


TIMES OY x 3 i) 


Key Rone xy XA 


bi ae ales 0° 3 


The table of values used in the tmnultiplication routine con- 
Siva ot a 6 x 16 matrix of product values ® Each product 
value is an exact value of a multiplication of two single 
hex-digit numbers. The SCaenae would normally contain the 
pmeoducus Of a zero multiplication 1S used for instructions 
within the ROM containing the table. A test for a zero input 
is made in the body of the subroutine. Each product value 
within the table consists of a two hex-digit number with 

the second value being the least significant digit. During 
the execution of the program, each digit in the two hex- 
digit number is loaded in a separate index register. In 
Tite mnette. wo. an M' “stands, form. the most “Significant digit, 
meio  Studanase or the least slenificant digit, and a "CY" 
Stands for a carry. The longhand multiplication with the 


DOreon Choppéedein this program can be shown as Lo]. lows: 
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3 3 1 0 
Yz Yop YY Y% 
ep ee 
m1 L3, M54 — —_— == 
Map bao Mop bap Mp — — 
M 
Mog Lg, M53 bog Mag M3 M3 — 
CY), CY. CY. CY, 
A reese eb 
4s bo a 26 


U1 
= 
iy = ss 


En) Pevelopment Of the trogram 
The method used in this procedure becomes very 
complex because of the large number of separate hex-digits 
involved and the small number of index registers available 
fometone eachwdadlett. This problem becomes more complex 
since each addition of two hex-digits creates a possible 
carrye & Solution to this problem is to make a process graph 


that simulates the multiplication process. 


i 





The process graph for the multiplication routine 


is shown in Figure 9. The symbols used have the following 


meaning: 
a Xo ~ table value of product of the 
digits X, times Y 
0 5 
hex-digit cy 
hex-digrit as sum addition 


The following examples is used to show the flow 


bf the hex-digits throughout the routine. 


Pieut. XxX = 
Y= 


Caleulavaon by) hand: 


EO eo 
cmon 
0 8 
Cf OnE 
Gem OC. is 2 
Cope a eee es 


cy cY cy CY 





Figure 10 shows the complete flow of the numbers throughout 
the routine. The value obtained agress with that computed 


by hand and can be used as a check of the process Pia. 
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SAMPLE MULTIPLY USING PROCESS GRAPH 


Figure \0o. 





Game COnsv~ructing the Program 

The process graph was used to assign and keep 
track of the Index Registers being used (Fig. 11). 

MicmsUbMOouUpINGcmrS: Stores all the values in 
the Index Registers into RAM, R2 and R3 are used to address 
values in the RAM. RO ana Ri are used to fetch values from 
the multiplication table. R4 > RF are the working registers 
and are assigned as shown in Figure ll. 

e. Expansion of Multiplication Routine 

The multiplication routine was written to give 
a truncated product of two positive numbers accurate to 
meer sipnaficant digits. The multiplication routine which 
handles both positive and negative numbers is shown in the 
process graph, Figure 12. 

Bach numbpersInpuvTed Into tne multiplication is 
meoved for its Sign. The sign bit is then shifted out to 
mermalize the number to four Significant digits. The pre- 
viously described multiplication routine is then executed. 
mae resulting number is shifted back to its proper form 
and the required sign bit is set. 

f. Summary 

The multiplication routine developed required 
the memory space of three ROMs and computed its result in 
merece mitliseconas. A compleue licting of the multiplication 


routine is found in Appendix C. 


at 
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MULTIPLY INDEX REGISTER ASSIGNMENT 


Figure \\. 
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2. The Cosine Routine 

The development of a Cosine routine for the micro- 
computer which could compute sufficiently fast was one of 
the major programming tasks. Cosine Routines written for 
general purpose computers are usually written as series 
approximations in order to save memory space. The object 
of the routine written for this program was to increase 
the speed of calculation. 

Two Cosine routines previously programmed on the 
MCS-4 Microcomputer were investigated. The first was a 
Chebyshev approximation routine which required 750 milli- 
seconds to compute the Cosine. The second routine investi- 
gated was a Cordic approximation which required 350 milli- 
seconds. Both methods were too time consuming for this 
project. 

The procedure developed in this project was a table 
look-up, linear interpolation routine. The Newton Divided- 
Difference interpolating Polynomial was used because of its 
Emote! ty. Lhe Size of the table réquired and the accuracy 
Stevunce resol! ts are DOUh/TunectvoOns Of the degree of the 
Polynomial used. For simplicity and speed, a first-order 
minmened—diiievence table was used which resulted in a linear 


itnverpolation of the form: 


F(X) = F(X) + (X - re, ee ae 
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where 


EG) = "COs. 06 


F(X, ) SO) 


Noe 


FLX, ,X_] = 
1"e i 0 


The table consisted of all values of F(X, ) and 
FLX 4, X54 for F(X) = COS 8, 0 < 6 < 1.88 radians in hexa- 
gdécimal. The increment used for each value in the table 
was .08 radians hexadecimal. The table was constructed 
from a Fortran program which used a decimal increment of 
Sse , equivalent to .08 hexadecimal, and outputted the 
desire table values in hexadecimal. 
2a ble Skene rh 
The size of the table loaded into the program 
was a function of the required accuracy of the Cosine rou- 
tine. Ihe data supplied to the program from the navigation 
Oeviceemuas cCccurate tO three significant figures. 
Diiterent size tables were constructed and 
EGsvedmelor A2cciuracy., ~olnce the interpolation was linear, 
Ucar ceouncirroOr OCeCurred "al the mic=point between each 


Cable value. 


JO1 









Maximum Error 


Linear Approx. 


The table size was adjusted until the maximum 
error Was WLtUhin three units in the fourth Significant 
figure, thus guaranteeing three significant figure accuracy. 

The table values were loaded sequentially into 
ROM number three. The first entry was the value of F(X) 
followed in order by F(X, 5X); P(X,), Nee? & P(X) FOX 44 X)- 
The remaining part of the ROM was used for the interpolating 
routine. 

Dee Wesecripeion Ol the Program 
The Cosine routine was designed to be called 


imomeouraslae Oi ncme ROM, ~The Input resisters RC, RD, RE, 


RF with RF the least Significant figure. 


Examples: 
ne RD RE RF 
Celts 520) would be loaded 0 2 5 1 
8 = 0.748 would be loaded 8 y 7 
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The description of the programs function is best shown by 


an example. Let the input be 6 = 0.74. Therefore 


RD RD RE RF 
0000 O100 Cia 0000 
0 4 rf 0 


The first object of the program was to find the 


required values (F(X, ) F(X X,)) in the table. The ROM 


altel e 
Cammbe thouent of es a 16 by 16 matrix of & bit words: Each 
Dyce therefore Will hold two significant figures of each 
table value. Since each table value has four significant 
figures and two table values were required for each compu- 
tation, only four bytes needed to be retrieved aeonntthe 
table. If the four required bytes for each computation are 
PMOMeMb Oh as sa unde, the Cable could be theught of as rows 
Sreunits WLth four units in each row, The increment in 
radian value therefore between each row is .20 radians hexa- 


Gecimal. The required row in the table is easily found by 


Gividing the inputted angle by two. 


—— = 03A0 (hexadecimal ) 
The required row in this example is row three. To find the 


reguiared unit adm row three, the remainder of the inputted 


angle is divided by four. 
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OOAO _ 
a 0028 


The required unit in this example is unit two. The final 
remainder of the inputted angle, 0008, is equal to the 
difference (X - X,). 

In summation, the location of FLX, 44%, 4 zB) 
the table was found by simple SEGUE by two and four 
performed on the input and the aifference. Divisions by 
powers of two can be obtained by shifting operations. 

After the desired table values and difference 
were found, the multiplication routine was called which 
multiplied (X - X.) times F(X, 44 2X4] and stored the result 
in index registers RC to RF. The final computation took 
the value in RC to RF and subtracted it from F(X, ) and stored 
the result in index register RC, RD, RE, and RF to be re- 
turned to the executive program. The step by Step procedure 
of the Cosine routine is shown in the Cosine process Seay 
acure 2. 3% 

ec. Expansion of Routine 

The table used in the Cosine Routine was also 
the table required by a Sine routine. It was noted that 
this routine coulda also be used to find the sin 9 for 
Oa Clg 010 aly, subtracting the input © from 90 (1.88 rads-6 


rads) and using the same routine since cos(90 - emo e 
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CT sin3szy 


hay oweooeOdd ANESOD 


a at eee vests 
ANISOO 


(x) a 


Z 


LOS 





qd. Summary 
The memory space required by the Cosine routine 
was one ROM plus the space taken up by the subroutines called 
Esmthne cosime routine. The main part of the Cosine routine 
contains only 46 instructions. The time required to execute 
the Cosine routine was basically the time required to execute 
the multiplication. The Cosine routine required only a 
total of 5.17 milliseconds. This computational Speed repre- 
sents a 70 fold decrease in the computational time to com- 
pute the Cosine by previously aval lanes roub@ne., spy cable 
look-up schemes, it was proven that the computational speed 
of the microcomputer could be competitive with that era 
general purpose computer. The fecic lpi las verot USI. tne 
mcS-4 as a navigation computer was also proven by Gemonstra- 
Piet seapeliny = boriindg the Cosine and Sine at competitive 
speeds and within the memory limits and accuracy limits 
required. 
3. Common Routines 

The common subroutines were written to do the basic 
housekeeping operations such as storing data, simple arith- 
metic, shift operations, and transfer of data between RAM 
and IR. These Subroutines were called by the executive 
routine and the multiply and Cosines routines to aid in the 
data handling. The functions handled by the common subrou- 
tines were those best suited for the MCS-4 and therefore 


could be written in a straight forward way requiring little 
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speed or memory space. The functions of the common routines 
were broken into three groups: Arithmetic, Shifting, and 
Data Handling. 

a. Arithmetic 

The arithmetic routines handled the simple addi- 
tions and subtractions required in the navigation program. 
These routines were handled well in the MCS-4 by the 4-bit 
ripple-through carry type adder incorporated in the 4OO4 
CPU. This allowed direct addition or subtraction of two 
hexydigits in either the accumulator Or RA [here were 
two addition routines, two subtraction routines, and two 
special purpose arithmetic routines written. 

The two addition routines written were ADDRAM 
and ADDRAMIR. ADDRAM adds the contents of the RAM addressed 
by IRO and IR1 to the contents of the RAM addressed by IR2 
and IR3 and stores the results in the RAM addressed by Live 
and IR3. The contents of the RAM addressed by TRO Vand Lal 


remains unchanged. 


RAM (IRO, IR1) RAM (IR2, IR3) 
RAM (IR2, IR3) 


ADDRAM occupies 12 words of ROM and takes 288.8 microseconds 





of computation time. 
ADDRAMIR adds the contents of the RAM addressed 


Due lnv@eddetiecO Une: CONLcnus Of Pao tn rumah ano: SlOres 
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Baee results tMmeinGe thru IRF. The contents of the RAM 


addressed by IRO and IRI] remains unchanged. 


RAM (IRO, IR1) 
! TR(C + F) 
IR(C > F) 


ADDRAMIR occupies 22 words of ROM and takes 237.6 micro- 
seconds of computation time. 

Thesuwo  SUubturactuon reutines written for the 
navigation program are SUBIR and SUBRAMIR. 

SUCTEReoUo tracer ote eomrenbs oil hG thru TRE 
from the contents of IR4 thru IR7 and stores the result in 


IRC thru IRF. The contents of IR4 thru IR7 remains unchanged. 


IR(4 > 7) 
TR(C + F) 
TR(C + F) 


SUBIR Occupies 17 words of ROM and takes 182.6 
microseconds of computational time. 

SUBRAMIR Mone naets ehescemvenus Of TRC thru IRE 
from the contents of the RAM addressed by IR2 and IR3 and 
stores the result in the RAM addressed by JIR2 and IR3. The 


contents of IRC thru IRF remains unchanged. 


RAM(IR2, IR3)———— 
RAM(IR2, IR3) 


IR(C > F) ———_-___> 
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SUBRAMIR occupies 32 bytes of ROM and takes 
345.6 microseconds of computational time. 

The two special arithmetic routines written for 
the navigation program are COMPLEMENT and COMANGLE. 

COMPLEMENT takes a four character hex number 
addressed by IRe and IR3 in RAM, takes the ones complement 
Gre eheo number, and stores the result back into the location 


in RAM addressed by IR2 and IR3. 


RAM(IR2, IR3) EN RAM(IR2, IR3) 


COMPLEMENT occupies 28 bytes of ROM and takes 302.4 micro- 
seconds of computational time. 

COMANGLE takes the angle loaded in IRC thru 
IRF, subtracts it from 90 in hex radians, and stores the 


mMeslvltepaeck-inbo IRC thru-lRr. 


TRC(C + F) TR(C > F) 


COMANGLE occupies 18 bytes of ROM and takes 194.4 micro- 
SecOmas Of computational time. 
Ben olen bite 
MilGiplvecation- and division of héey=difrits by a 
multiple of two was accomplished by shifting the variable 


Cicer Ct trom Pach: thiesreaquired number of bits. To take 
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advantage of this capability, two subroutines were written 
HOimEne Navagation program, DIV2iR and D1Ve2. 

DivV2ZIR divides the contents of IRC thru IRF by 
two and stores the result in the RAM addressed by IRO and 
IR1. The division iS accomplished by shifting the variable 
Peeno One bee, Ihe original contents of TRC thru IRF are 


unaffected. 


IR(C + F) Lo RAM(IRO, IR1) 
Th 


DIV2IR occupies 24 bytes of ROM and takes 259.2 microseconds 
of computational time. 

DIiVe divides the contents of the RAM addressed 
by IR2 and IR3 by two and stores the result in the RAM 


addressed by IR2 and IR3. 


RAM(IR2, IR3) a RAM(IR2, IR3) 


DIV2 occupies 33 bytes of ROM and takes 356.4 microseconds 
of computational time. 
Cae Data Handling 
The transfer or data within the navigation pro- 
gram was accomplished by the data handling routines. These 
routines saved much memory space in the Executive Routine 


and the Multiply and Cosine Routines by grouping, these tasks 
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into separate subroutine calls. There were five data handling 
routines written for the Navigation program. Three routines, 
RAMIRC, RAMIR8, and IRRAMC, were written to transfer data 
between the IRs and RAM. TRANRAM was written to transfer 
data between different locations in RAM. TIME was a special 
routine written to load the proper time interval of one navi- 
Papronecycle’ inco the IRs. All five routines were written 
fo handle data of four hex=-digit size. 

RAMIRC transfers the contents of ane RAM addressed 
Praite song sths into IRC thru IR. The original contents 


of RAM are unaffected. 


RAM(IR2, IR3) — IR(C + F) 


RAMIRC occupies 16 bytes of ROM and takes 172.8 microseconds 
SusconpuLtational time. 

RAMIR8 transfers the contents of the RAM addressed 
Dy Uno and UR into IRS thru IRB: The original contents 


of RAM are unaffected. 
RAM(IRO, IR1) —> (IR(8 + B) 


RAMIR8 occupies 16 bytes of ROM and takes 172.8 microseconds 


of computational time. 


amas 





IRRAMC transfers the contents of the IRC thru 
Piveencremricennwmeaddressed by ER2 and IR3. The original 


contents of IRC thru IRF are unaffected. 


IEC ae" le), —> peavic ihe sls) 


TRRAMC occupies 16 bytes of ROM and takes 172.8 microseconds 
of computational time. 

Tro habineroansters ac, COmuemus Ol  hAM addressed 
by IRe and IR3 into RAM addressed by IRO and IR1. The 
original contents of RAM addressed by IR2 and IR3 are 


unaffected. 


RAM(IR2, IR3) sos) RAM(IRO,IR1) 


TRANRAM occupies 8 bytes of ROM and takes 86.4 microseconds 
of computational time. 

tiie see special suproutane wee Lven to load 
the time interval of one navigation cycle into IR8 thru IRB. 
The time interval used in this program was 200 milliseconds. 
Vominpus a diiferent Gime ianterval into the navigation pro- 


gram, a nev TIME routine is substituted into the program. 


—— + B) 


dbdee 





TIME occupies 9 bytes of ROM and takes 97.2 microseconds 
of computational time. 

ec. Summary 

Thirteen Common Subroutines were written for 
the Navigation program. Each routine was designed to 
handle the particular needs of the Navigation program. 
All *enerveen routines fit in one ROM. A complete listing 
of the Common Subroutines is in Appendix A. 
4. Summary of Subroutines 

Each operation defined in the operational process 
graph, Figure 6, was successfully programmed within the 
memory and time constraint of the program analysis. Each 
PPL OubIne Was written in a modular form allowing for easy 
aedicgion Or subtraction of new code. The memory size, 
computational time, and capability of each subroutine 1s 


impsted in table Vil. 


Pe ecUrIVvE ROUTINE 

The Executive Routine was written to call up the sub- 
routines in the order described by the operational process 
graph, Figure 6. The Executive Routine established the 
priorities of each function and was designed to make all 
the decisions in the execution of the Navigation program. 
The variables used by the Executive Routine were all Le onea 
in RAM. The designated location of each variable in RAM 


MS chown on Taole VIII. 


dels 





Subroutine 


MULT 

COS 
ADDRAM 
ADDRAMIR 
peel 
SUBRAMIR 
COMPLEMENT 
COMANGLE 
DIV2IR 
ede 
RAMIRC 
RAMIR& 
IRRAMC 
TRANRAM 


TIME 


Length 
(8-bit 
words ) 


743 
ae 
12 
22 
17 


32 


18 
Ou 
oo 
ne 
16 


16 


Table VII 


Time 


(usec) 


5000 


E70 


289 


238 


183 


346 


502 


194 


25) 


356 


iS 


a 


ii 


86 


JI 


Purpose 


[IR(8+B) ]- [IR(C+F) J*+IR(C+F) 
COS[IR(C+F)] + IR(C+F) 
RAM(0,1) + RAM(2,3)+RAM(2, 3) 
RAM(0,1) + IR(C+F)+IR(C+F) 
ot Ro 
RAM(2,3)-IR(C+F)+RAM(2, 3) 
Sat Tec may) alee 1 3)) 

(90° )-IR(C-F)+IR(C+F) 

PIRK CE) 172>RAM(O, 1) 
[RAM(2,3) ]/2>RAM(2, 3) 
RAM(2,3) + IR(C+F) 

RAM(0,1) + IR(8+B) 

IR(C+F) + RAM(2,3) 

RAM(2,3) + RAM(0,1) 


(Nav Cycle Time)~>IR(67B) 


NAVICATION SUBROUTINES 
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Mme iansy Operation of the Executive Routine was to 
compute the proper Cosine and Sine of the current true 
heading. Before calling the COSINE subroutine, the Execu- 
tive Routine determined which quadrant the true heading was 
in and set up the input angle in the proper form so that the 
Sine or Cosine could be computed by the COSINE subroutine. 
After calling the COSINE subroutine, the Executive Routine 
SeOoread the Output with the proper Sign in its location in 
RAM. A graphical display of the decisions made by the Execu- 
tive Routine in computing the Cosine and Sine of the true 
heading is illustrated in Figure 14. 

Next the Executive Routine Smoothed the current TAS 
halt Catlony With the previous TAS computed. The smoothed 
TAS was then used to compute VAX and VAY as described in 
EMiemNaviteatlcn equabions and illustrated by the operational 
process graph. 

After computing VAX and VAY, the Executive Routine tested 
to see if the INS was good. If the INS was good, the Execu- 
tive Routine computed the wind acting on the vehicle and 
the distance increments traveled as measured by the INS. 

If the INS was down, the Executive Routine computed the dis- 
tance increments traveled as measured by the AIR-MASS system 
using the last computed wind components. Next the Doppler 
System was tested to see if it was up or down. If the 
Pepper was up, the distance increments as measured by 


the Doppler were computed. 
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The distance increments outputted from the Navigation 
Program were selected by the Executive Routine in the 


Polovwene Ppriormey according to the following conditions: 


Pieisogei ty Conditions 
ite tNS and Doppler INS Upsamd Doppler Up 
Smoothed 
Comlidl > INS up and Doppler down 
Baan DODPLEr INS down and Doppler up 
4, Air-Mass INS down and Doppler down 


The Executive Routine was loaded into two ROMs with 
space left for addition of new code. A complete listing 


Of the Exeuctive Routine is in Appendix B. 


PY. ERROR ANALYSIS 

Microcomputers have a limited arithmetic capability. 
The basic arithmetic operation of the MCS-4 is a four bit 
addition. More complex arithmetic processes had to be 
reduced to this basic operation. It was therefore very 
MnmoOrPEAMNEL CO avold Unnecessary precision throughout the 
SalemEotloOnoesotnce the s2npubs anto. the Navigation program 
came from instruments whose precision is limited to three 
hexadecimal digits, the choice of four hexadecimal arith- 
metic waS conSidered to be Sufficiently accurate. The 
following, error bound analysis was performed to show that 
the input errors dominate the total error. 

The starting point for the error analysis is the opera- 


tional process praph, Figure 6G. It was apparent from the 
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operational process graph that the outputs DX and DY are 
symmetric, therefore an analysis of only the computations 
for DX were made. A process graph that involved only the 
operations which have an influence on DX was constructed, 
as shown in Figure 15, from the operational process graph, 
Figure 6. 

‘The errors corresponding to each operation in Figure 
were designated CE ,> €o> C35 seey Ep. The initial errors of 
the inputted data were expressed as e(TH), e(VGXI), e(SDD), 
and e(SHD). Due to the smallness of the errors, the products 
of errors were considered negligible when compared to the 
linear terms. The propagation error was derived as discussed 
in Chapter 2 of Reference 28 in the Bibliography. 

There were two means by which each operation contributed 
fo the error propagation. 


1.) Transmitting the errors which were inputted into 
the operation. 


2.) Adding an error of its own, which is due to the 
rounding or truncating process which limits the 
number of digits carried to the next operation. 


The transmitted errors were calculated by calculating 


Themotiterengials «oL the expression. 


dirty) = sax + dy 


disc m=) dt x dy 


X ydax-—- x dy 
Age eee 
VY Cc 
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a(sin x) = cos x ax 


a(cos x) ~ sin x ax 


The rounding or truncating errors were simply added to 
the transmitted error and thereafter propagated through the 
Pemedinder of the calculation. 

The error after the first operation in Figure 15 is 


given by: 


Secos Cll) ia Sel CLM ene ey 


This error is further transmitted bv operation three. 


e(cos(TH)SDD) = SDD(-sin(TH)e(TH) + e,) + cos(TH)e(SDD) 
+ 


oe 


pimilarly, the error alter the second and fourth operation 


is @iven by: 


e(sin(TH)SHD) = SHD(cos(TH)e(TH) + e.) + sin(TH)e(SHD) 


a ey 


Prem grvess thneetotal error alter OpGieayi0on Six as: 


te 





e(cos(TH)SDD) + e(sin(TH)SHD) + e¢ 
In a similar fashion, the error after operation five is: 
e(T-VGXI) = VGXI-e(T) + T-e(VGXI) + C5 
Combining the above, the total error is expressed as: 


e (DX) = ={e(T-VGXI) + [e(cos(TH)SDD)+e(sin(TH) )SHD + eg te tte, 


The following is obtained by expanding the expression for 


e(Dx): 


e(Dx) s{VGXI-+e(T)+e(VGXI)- Tre +SDD(~sin(TH)e(TH)+e, ) 


D 


+ cos(TH)-e(SDD)te +SHD(cos(TH)e(TH)te,) 


3 


+ sin(TH)e(SHD)+e,+tegte.teg 


The error bound for e(Dx) is computed by using the triangular 


inequalities. 


|e (Dx) | 


lA 


S{|VGxI|- er) | + |T|-|e(VGxI)| + |SDD|[]sin(TH)| (TH) 


+ |e, |] + |cos(TH)| |e(SpD) [+ | SHD| [|] cos(TH)| |e CTH) | 
= leone) a enn CDH ie (SHD) || > + Jec| + Je. Penn 
+ 


leg] + lel} + leg 


eae 





The time increment T is controlled by the programmer 
and hence may be chosen to be a power of sixteen, therefore 
the inherent error e(T) = 0. The addition process carried 
out in operations six and seven is carried out exactly in 
the range of values permitted, therefore Ee 0 and en = Oe 


This reduces the error bound expression for the total 


Preperaved eCLreie CO; 


e(Dx) < 5{|T|-|e(vexr) | + |SpDD|[]sin(TH)|]e(TH)| + le, | 1 
+ |cos(TH)||e(SDD)|+ |SHD|[|cos(TH)||e(TH)| + [e,]1] 
+ peer letSHDy Ts lec | + le. mary ees leg | 


ine error pound for each operation in Figure 15 due Go 


roundoff and truncation was found to be: 


Operation Error Bound (Decimal) Corres. Error 
Cosine a 1072 ey 

Sine Bax 10_¢ C5 
Multiplication See Oy e5 fey ls lect 
Payals ion 8 x 10_¢ ep 
Subtraction 8 x 10 

Addition 0 |e 6 : eo 








Theverror DOUnGa 1Or tne inputs was Obtained from published 
Seu ces wero imceran ACLUG es Mmaximumecerror Tor -each system could 


Nev be foun sene (30) value, 99.7e CEP, was used. The 


lee s5 





systems onboard the P3C naval aircraft were used as repre- 
senting the current "state-of-the-art" systems in operational 


use today. 


System Error (30) Corresponding Error 
Inertial Navigation 5.4 NM/HR |e (VGXT) | 
System 
Inertial True Heading .5 degrees |e (TH) | 
Doppler Along Heading 1.8 |e(SHD) | 
Doppler Across Heading 3.6 |e (SDD) | 


The error bound for the Doppler is a function of velocity. 
ineorder to compute the error bound for the system, an exam— 
ple was used. The error bound analysis was calculated for 
pmoereratt) with a velocity of 400 knots along the true 
heading and a velocity of 50 knots across the true heading. 
The navigation cycle time used was 150 milliseconds. The 
true heading of the vehicle was 045 degrees. The inputs 
and error bounds inputted into the system were calculated 


as follows: 


TH = O45 degrees = 0.7854 radians 


e(TH) = .5 degrees = 0.0087266 


VGXI = 3.8.2 knots 


e(VGXI) = 5.4 knots 


Tite = .15 sec = .0000416 ER 
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SHD = (400 knots)(.0000416 HR) = .01664 mm. 
2 (cotDee=9 4. Oloo" nm)y(,010) =). 0002995 nm. 
SDD = (50 knots)(.0000416 HR) = .00208 nm. 


e(SDD) = (.00208 nm)(.036) = .0000748 nm. 


Although the values calculated above are dimensionally 
eorrect, they are not in the proper form for inserting into 
the onOueeT. The MCS-4 microcomputer used in this system 
works with a 16-bit fixed point variable. In order not to 
loser sitenitftcans Gigits in each number, the inputs are all 
normalized. The fixed point number used in the calculation 
Mom demrOrmeon 1°), bierefore the value of VGXI 
used in the computations is .3182. The time variable is 
normalized to .4160 which results in SHD equalling .1664. 
Now that the inputs are dimensionally correct and in the 
Beoper toOrmatv, they can be Substituted into the error bound 


equation as follows: 
fe:@Des)) |) 5{(.416)(.0054) + (.0208)[(SIN(.7854))(.0087266) 
TO OeO 2) (COs (7854) ) (2000746) 4+ (.1664) 
CC Oe 7054) (20067205) + 200003] + (SIN(.7854)) 


(.002995) + .00003 + .00003 + .00003} + .000008 
le(Dx)] < .0030798 





This is the maximum error bound for the total naviga- 
tion system. The maximum error bound equation was broken 
into two parts, one part indicated the maximum propagated 
errors due to the inputs and the other indicated the maximum 
error developed by the navigation program. 

The error bound for the inputs alone is ‘calculated as 


MOMmLOWS : 


le(Dx) Ts put < 5{|T| |e (vex) | 7 SDD i Sane Guns recrH)) | 


+ \COS (TH) | le(Spp)| + |SHmi|CcOS (TH) |e(TH) 1 


Pet CEH) ieccHy) |p 


2 - 
le (D4) Tan put < 5{(.416)(.0054) + (.208)(.70711) (. 0087266) 


revOy i) (eOO0R4 ep Ge MOG Gry Oro ds) C0007 20eD 


a 


te 7O7a) (2002095) 3 


Je (Dx) put = ,003024 


Nest the maximum error due to the calculations in the 


ievawerion program alone were. computed. 
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le(Dx)| < ${(.0208)(.00003) + (.1664)(.00003) + .00003 


prog — 


+ ,00003 + .00003} + .000008 


[et2b sy eyes < .0000557 
The results of the error bound analysis indicated that 
the maximum error created in the navigation microcomputers 
computations was only 1.8 per cent of the total maximum 
error. It was concluded from this analysis that the accu- 
racy of the MCS-4 navigation program using a 16-bit fixed 
Meraeadtua Lenmeth wes well) within the limits required for 
the navigation problem. 

The best way to check the results of this error bound 
analysis would be to fly the system in an actual aircraft. 
Sameer an alrcratt was not available for this purpose, a 
detailed FORTRAN simulation program was written to test 
Preminimetitons of the navigation program. Section VI dis- 
cusses the FORTRAN simulation program and the results 


obtained from examples tested on this program. 


G. SUMMARY OF PROGRAM 


The navigation program written in this thesis for the 


MCS-4 microcomputer takes the outputs of the Inertial navi- 


gation system, Doppler navigation systen, and Air-Data 


he 





system and computes the vehicle's distance traveled ial ela 
east/west and north/south directions. The total program 
consists of 1768 instruction words on seven 4001 ROM chips 
broken into an executive routine and fifteen subroutines. 
Two 4002 RAM chips are required to store the data and varia- 
bles used in this program. The total computational time 
required for one navigational cycle is between 36 and 80 
milliseconds depending on the navigational mode used. The 
computational error developed by the navigation program from 
error bound analysis represents only .1 per cent of the 
total error. The program uses a sixteen bit fixed Dol mer 
yariable which allows it to accept inputs up to 2047.99 
knots with an accuracy of + .0625 kts. The total cost of 
the one CPU chip, seven ROMs, and two RAMs used by the system 
ve 095,00. 

A complete listing of the navigation program written 
in the MCS-4 assembler language is found in Appendices A, 


Beeana C. 
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VI. FORTRAN SIMULATION PROGRAM 


The development of the MCS-4 Microcomputer navigation 
program was a very time-consuming process that involved the 
construction and testing of many subroutines. The testing 
of the total program in its final form was jameced by. eine 
capability of modeling the inputs to the navigation program 
in a machine level language. Hand calculated examples were 
tested and used to debug the completed program, it was de- 
cided to write a detailed bit-by-—bit digital simulation of 
the navigation program in FORTRAN on the IBM 360/67 compu- 
ter. The FORTRAN simulation was written so that major 
changes in word length and scaling could be tested without 
requiring major rewriting of the microcomputer program. 

Parameter changes were made to the inert ual poppilews 
and Air-Mass inputs without incurring the added cost and 
time of developing an actual system. The FORTRAN simulation 
was a comprehensive and flexible means OPrTtesting swe Macres 
computer program through a series Of realistic and unrealis— 
tic tasks. A complete listing of the FORTRAN simulation 


program is found in AppenGaw). 


A. DEVELOPMENT 
The FORTRAN simulation program was developed aS an exact 
simulation of the navigation program developed for the MCS-4 


microcomputer. The program was written by utilizing the 


process graphs developed in Section V. The program includes 
a parallel solution of the navigation equations utilizing 
the FORTRAN routines available on the IBM 360/67 computer. 
fhe results Of the solution by the MCS—4 simulation and the 
FORTRAN solution provided a comparison of: the 16-bit fixed 
point solution of the MCS-4 microcomputer program and the 
24-pit floating point routines on the IBM general purpose 
CoOmoucer. 

The FORTRAN Simulation provided the programmer with the 
ability to input many variables. The outputs of the Iner- 
tial, Doppler, and Air-Mass systems were modeled to include 
system failures and input errors. The outputs of the simu- 
lation indicated how well the programmed microcomputer could 
handle changing variables as compared to a programmed 
general purpose computer. 

The main body of the simulation program was written as 
a subroutine in order that systems programmed or simulated 
in FORTRAN could easily be tested with the navigation routine. 
The main body of the simulation program declares the index 
registers, random access memory locations, and read only 
memory locations as dimensioned arrays. The numbering and 
size of each array corresponds exactly with the addressing 
and size of these areas within the MCS-4. The main body of 
the simulation program prepares the inputs and initial con- 
ditions for the execution of a simulated navigation cycle. 


Piewiera@eeimal Cosine table stored In the MCcS-!| ROM is 
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computed and stored in the ROM array by the main body of 
the simulation program. The simulation begins when the 
Subroutine NAV is called with the initial inputs. 

Subroutine NAV is a FORTRAN simulation of the MCS-—4 
navigation executive program. The inputted variables to 
NAV are 32-bit floating point decimal numbers. Each varia- 
ble is converted to a 16 bit fixed point hexadecimal number 
By the subroutine CONVRT, Appendix D. After each input is 
in the proper form, subroutine NAV stores each variable in 
Belleman array at the same location corresponding tO the RAM 
addresses in the MCS-4, Table VIII. After the inputs have 
been stored in the proper locations and in the proper format, 
subroutine NAV ealls a FORTRAN simulation of each subroutine 
called by the MCS-4 navigation executive routine. Each 
subroutine is executed in the same order as in the MCS-4 
executive program. The 16-bit hexadecimal results are 
printed to provide a check with the results of the actual 
MCS-4 microcomputer navigation program. The NAV subroutine 
iaeriyvecalls subroutine CONVD, Appendix DP, which converts 
the 16 bit hexadecimal results into floating point decimal 
results which are then compared with the solution of the 
navigation equations by floating point FORTRAN calculations. 
Ppeonobete Listing of subroutine NAY 1s found in Appendix D. 

Each common subroutine programmed on the MCS-4 micro- 
computer was duplicated in the FORTRAN simulation. A conm- 


Dlete description of each subroutine and ite graphical 


tec el 





representation can be found in Section V.D.3. A complete 
listing of the FORTRAN simulation of each of the common 
Subroutines is found in Appendix D. 

The multiplication and Cosine subroutines used in the 
FORTRAN simulation were developed from the multiplication 
process graph, Figure 9, and the Cosine Process praph., 
Figure 13. The Cosine routine uses the Cosine table stored 
by the main body of the simulation program. The multiplica- 
tion table used by the simulated multiply routine was com- 
muted individually for each pair of hexadecimal digits 
pailed by the multiply routine. 

The last subroutine written for the simulation program, 
subroutine PRINT, was written to incorporate some of the 
features of the MCS-4 Assembler and Interpreter. Subroutine 
PRINT, Appendix D, was developed to allow the programmer to 
dump the contents of the index registers and RAM locations 
at any point during the execution of the simulation. 

A complete listing of the FORTRAN simulation program 


is in Appendix D. 


B. APPLICATION 

The FORTRAN simulation was first used to optimize the 
navigation program. The first favipotion program Wratten 
on the MCS-!| had been designed to handle only positive 
numbers within the true heading limitations of zero to ninety 


degrees. Before expanding the propram to handle all quadrants 





and both positive and negative inputs, different ways to 
implement these changes were investigated. The Fortran 
simulation allowed numerous techniques to be attempted with 
the advantages and limitations of each readily apparent 

in the output. The Fortran simulation program was easy tO 
change, whereas the MCS-4 program would have required major 
program changes to implement each @echnique. After the 
Fortran program had been modified to do the required tasks, 
the MCS-4 navigation program was modified to incorporate 
these required changes. The true heading was used in four 
separate quadrants of ninety degrees each with the Execu- 
tive routine testing the true heading input and computing 
the Cosine and Sine as required. The sign of each variable 
was incorporated into the program by using a ones complement 
scheme with the first bit of the 16-bit word becoming the 
sign bit. The required changes to each subroutine was first 
tested in the simulation program before being programmed 
jmee the navigation routine. 

The second use of the simulation program was to debug 
the MCS-4 program in its final form. The Simulation program 
was run through a series of examples. The simulated MCS-4 
program results were compared Valais tne POruran compured 
results to insure the correctness of the Simulation program. 
After the simulation program had been completely debugged, 
the actual MCS-4 program was run with the same inputs as the 


simulation prorram. By dumping the contents of the registers 





and RAM memory at key points in the actual program and 
comparing these values with those computed and printed out 
by the Simulation program, the sections of the actual pro- 
gram in error were easily identified. After changing a 
section in error, the actual program and simulation pro- 
gram were run to see if the correction was correct and also 
see if it affected any of the other sections. This method 
of checking and rechecking proved to be an excellent tool 
for debugging microcomputer programs. The only problem 
encountered in this method was the required time for outputs 
to be printed from the IBM 360/67. 

Recmeneira wse Ol the Simulation program was to check 
the results of the error bound analysis. The same inputs 
moea in the error analysis Were inputted into the simulation. 
Simulation runs were made uSing correct invduts in order to 
compare the results of the MCS-4 calculation with those of 
fe FOrtran program in order to determine the errors due to 
ice navigation program. Simulation runs were made for each 
Mavigation mode of operation and for each true heading quad- 
rant. Next the same runs were made with the inputs at 
maximum errors. The results of these runs and the comparison 
of the errors developed are summarized in Table IX. 

The results, summarized in Table IX, confirmed the re- 
Sults of the error bound analysis. It was noted that the 
ereawest computational ¢rror occurred when the vehicle tra- 


veled a direct path with constant inputs. This was due to 
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the linear addition of the truncation error when the inputs 
remain constant. It was also noted that the computational 
error was zero when the vehicle returned to the departing 

position indicating that the truncation error cancelled in 


the opposite direction. 
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VII. CONCLUSIONS 


The results of this design study indicated that a micro- 
computer is both fast enough and powerful enough to handle 
the complex task of navigation. The total computational 
time of the navigational cycle, 80 milliseconds, is well 
below the 200 millisecond cycle time uSed in current naviga- 
tion computers. Table look-up routines for such complex 
tasks as multiply, Cosine, and Sine proved an effective 
means of trading inexpensive memory for computational Speed. 
The total hardware cost of the ten LSI cips, excluding a 
board to hold the DIP packages, was $95.00. 

The software aids investigated in this study were very 
effective in decreasing programming and debugging time. 
Ceapmlneory, in the form of process graphs, was @n excellent 
means of visualizing the actual flow of the data throughout 
the computations. The initial problem was succesSively 
broken into smaller and smaller discrete parts until each 
operation could be programmed easily in the assembly language 
of the MCS-4. The operational process graph was then used 
£0 combine the proprammed discrete operations into the re- 
quired navigation computations. The actual writing of the 
program was greatly simplified by this method. 

The error bound analysis of the microprogram was greatly 
simplified by the uSe of proces& fraphs. The maximum error 


jrmeeren of the diserete operation&® Was initially computed. 
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The summing of these errors as they propagated throughout 

the program was easily identified with the aid of the opera- 
tional process graph. The testing of actual erros in the 
computation, however, proved very difficult. Sample problems 
were hand calculated with the aid of the process graphs and 
the results were compared with the results computed by the 
actual program. When these results did not agree, both the 
hand calculated procedure and the computer program had to 

be debugged. This procedure was very time consuming with 
only a limited number of tests being made. The FORTRAN 
Simulation program was written to aid in testing and debugging 
the navigation program. 

The FORTRAN simulation program proved to be an effective 
erred LeStane the navigation program. The set up of the 
srmubataon program allowed the results of the navigation 
program to be compared at different points in the computa- 
tions with a FORTRAN solution of the same required equations. 
With the aid of the IBM 360/67 computer, many tests were run 
on the FORTRAN simulation program. The results of these 
runs identified those areas in the navigation program that 
needed rewriting. Before expanding the navigation program 
to include negative inputs from different true heading quad- 
rants, the FORTRAN simulation program was changed to incor- 
porate the required changes. After all the changes had been 
tested and optimized on the simulation program, the MCS-4 


Navigation program was rewritten to incorporate those changes. 





The simulation program proved very effective in optimizing, 
expanding, debugging, and testing the navigation program. 

It is concluded from this design study that many of 
the dedicated computational tasks being done by large 
general purpose computers on board naval ships and aircraft 
can be done by microcomputers. It is recommended that 
research and development of a multi-microcomputer system be 
begun to replace one or more of the costly general purpose 
computer systems currently being used. The development of 
the MCS-8080 microcomputer with its increased computational 
speed and power over the MCS-4 will make this task easier. 

It is also recommended that before programming a@ micro- 
computer, the desired program should be written in a higher 
level language simulation program before investing in the 
software costs of a microprogram. The development of PL/M, 
ee from PL/I, a high level language for the MCS-8 ana 
MCS-8080 greatly increases the ease of programming and 
testing programs for the microcomputer. 

The rapid improvements being made in microcomputers, 


and the current availability of inexpensive microcomputers, 


make it imperative that the Navy begin now in the development 


of systems utilizing this technology. 
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Cie. 
SRG R2$ 
RDM; 
CMA; 
TAC $ 
WRI4 3 
INC R33 
SRC P23 
EDM 20 : 
XCH ROQ3 
RDM; 
CMA3 
ADD RQ; 
WRIM5 
INC R33 
BG Re 5 
RDIM3 
CMA; 
ADD RQ; 
WRMs$ 
INC R33 
SEG R23 
RDM; 
CHAS 
ADB &Q3 
WRM 35 
CLB; 
Sow 70; 
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CMC; 
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CHC; 
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XCH RF; 
CLB; 
BBL QO; | 
ENDSUB: NOP; 
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APPENDIX B 
Me AEN He A BE Ae 2 oh a Xe 2 Shea 2 aK ee aes ae a a ae a eo a a a ee ak aft ea Kk ake ak ake a Ae kak ok a ke ae 2k ake ok 
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APPENDIX ¢ 


she oh a sig a ak aft ok ake ake ae ak ake ok ake ake Se ae ak a ae ke a ate ae oe te fe ak ko ae oe ate abe te a aR fe at a oe aie abe ob ak ake ok ake ak ate ok 
THIS SECTION CONTAINS THE NAVIGATION SUBRCUTINES 
af se 2k i fe oi 2K eo a oe a eo ae ec a oe he oe a ok ok ake ak oe ac ak oa ke oe oe ak aK ke ake ok Xe oe ak cake eo ok 
BE SR TR IK RE Be AR AE NE Ns AE Me Oe a EE ce ee ae ete ea oie Eo EE OE KO OK ee a OK ie oe oe ok OK 
THIS 1S THE TABLE OF VALUES USED BY THE COSINE ROUTINE 
THE FIRST COLUMN CONSISTS OF THE VALUE OF THE COSINE 
FROM 0 TO 90 DEGREES IN INCREMENTS OF .08 HEX RADIANS} 
THE SECOND COLUMN CONSISTS OF VALUES USED TO 
INTERPGLATE BETWEEN VALUES IN THE FIRST COLUMN; 
MO FE RE EAE AE AK Bs As RPS CE eC CAS AE EO CE AC AE OK AEE AE 2 oe EE OK KE SE OK OK eRe Me EE A EEE OE A eK 

CON "FEEE"; CON "0040"; 

CON “OEFE"; CON "FESO"; 

CON “OSFE"; CON "BF31"3 

CON "OEEF"; CON "2F81"3 

CON "LOEF"; CON "LE32"3 

CON “2ECEN; CON "8C32"; 

CON "38BE"; CON "4A33"3 

CON "GEESE"; CON "3782"; 

CON "BO8F"; CON "3334"; 

CON "1LESEN3 CON "3EA4";3 

CON “AQ3E"; CON "F725" 

CON "6Q1E"; CON "70A5"3 

CON "63EE"; CON "9716"; 

CON "A2BE"; CON "2D86"3 

CON “3E7E"5 CON "1107"; 

CCH "364"; CON "437743 

CON "SAOE"; CON WOZET7US 

CON "76CD"; CON "F158"; 

CGN “FRED'S CON "3EB8"; 

CON "F24D"3 CON "4829" 3 

CON “BOFC"; CON "1099; 

CON "3DAC"S CON “E5ES"3 

CON "9D5C"3 CON "685A"; 

CON "CAOC"S CON "C8BA"; 

CON "0582"; CON "6618"; 

CON "5C5B"; CON "417B"3 

CON "COOBN; CON "4SCB"5 

CON "82AAN; CON "SE1LC" 

CON "814A"; CON "607C"; 

CON "OEDO"; CON "SFBC"; 

COW "1879"; CON "1500" 

CON "BECO" ; CON "6350"; 

CON "15A8"; CON "EBSD"; 

CON "5828"; CON "7ADD"3 

CON "BSC7"3 CON "CELE; 

CON "S857"; CON "§35E"; 

CON "26E4"3 CON "BABE"; 

CON "C176" CON "GESEN; 

CON "DFS"; CON "2EEE"; 

CON "6485"; CON "GALE" 

CGN "S805"; CON "G24E" 

CON NE LO4G"S CON eo r 

CON "4614"; CON "D78F"; 

CUN "0A93"; CON "AGAENS 

CON "CC 13"; CON "ODRENS 

CON "GFS2"; CON "S20F"3 

CON "6022" CON "83EF"; 

CON "4LAL"; CON "9OFF'; 

CON wel 73 CON UAQEE Ms 

CON "FLAO"; CCN “BEFE"; 
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Be AK AS AR MK RE AR TE SS 2 Ae oe he Fee HE BK Ae BC a Ae Me a CAE AE a ae Re Oe he oie Me AE oie a ee As hE OK ES aS AE Oe EE oe 
THIS ROM CONTAINS A 16%16 HEX MULTIPLICATION TABLE. 
ESCnmvetuG REPRESENTS A MULTIPLICATICN OF 7WO SINGLE 
eee lt NUMBERS «THE ZERO ROW 15 USED FOR TABLE 
INSTRUCTIONS; 
He SE Be Ree RET AERC HE SEE RC RCE BO AR RC OR AE OK A AR OK ae ae a ae ak he ake atc ak ak ake ae ak a ae ak ae ai 2 ak ak leak 2 
TABLE: FIN RQ; BBL QO; NOP ; NOP; 
NOP 5 NOP ; NOP $ NOP ; 
NOP; NOP; NOP $ NOP 5 
NOP ; MOP 5 NOP; OP ; 
miuWesmecones OO > CON “OL, CEN “az, CCN “o3"s 
Conm oo; COUN MOa", iCON Olio"; “CGN "0 Ct. 
CoM OG CON ON, CON “OATS CON "OBNS 
Clee OC! CON: ODM. CON “OEM; CON “OE; 
menaseceN "OO"; CON “O02, CON "G45 CON "O06"; 
COs OSes GCUNnOn s CON UOC, CON POEs. 
CONmEO"; CONSE IT2"; CON “E4'"'s CON aaho",; 
PUN LOO Pas CON TC ss CON Tre, 
we. CONTE OOM ECON “O3". “CON "O06" 5 CON "“O9N; 
SUNG CON MOR CON Wh Z se CON halts 
Sel meow CON boll CON VLE", “CON M2its 
OMe Gi (CONVO 2 CON NZ AY CON Tea. 
icemeeenN OOS CON “Oats CON "O8N, CEN TOC. , 
CUNpe LOC SCON Pat CON ee CON re: 
CRN ZO" eC ON 24s CON Zo Se ON se 2c, 
CoN Ss Oe BOON 934 CON oo CUNT Ce. 
SUNS cone OO; CONS "O5", CON “OA aC Ole OF. 
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